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ABSTRACT 

We perform an extensive test of theoretical stellar models for main-sequence stars in ugriz, using cluster 
fiducial sequences obtained in the previous paper of this series. We generate a set of isochrones using the 
Yale Rotating Evolutionary Code (YREC) with updated input physics, and derive magnitudes and colors in 
ugriz from MARCS model atmospheres. These models match cluster main sequences over a wide range of 
metallicity within the errors of the adopted cluster parameters. However, we find a large discrepancy of model 
colors at the lower main sequence {T^a < 4500 K) for clusters at and above solar metallicity. We also reach 
similar conclusions using the theoretical isochrones of Girardi et al. and Dotter et al., but our new models are 
generally in better agreement with the data. Using our theoretical isochrones, we also derive main-sequence 
fitting distances and turn-off ages for five key globular clusters, and demonstrate the ability to derive these 
quantities from photometric data in the Sloan Digital Sky Survey. In particular, we exploit multiple color 
indices {g- r, g- i, and g-z) in the parameter estimation, which allows us to evaluate internal systematic 
errors. Our distance estimates, with an error of cr^m-M) = 0.03-0.1 1 mag for individual clusters, are consistent 
with Hipparcos-hased subdwarf fitting distances derived in the Johnson-Cousins or Stromgren photometric 
systems. 

Subject headings: globular clusters: individual (M3, M5, M13, M15, M71, M92) — Hertzsprung-Russell 
diagram — open clusters and associations: individual (M67, NGC 6791) — stars: evolution 
— Surveys 



1. INTRODUCTION 

Improved stellar distance estimates are required to address 
many outstanding problems concerning the formation and 
evolution of the Galaxy. However, reliable distances from 
trigono metric parallaxes, such as those from the Hipparcos 
mission (lESAl 119 97). are restricted to a few hundred parsecs 
from the Sun. Therefore, for now, we must rely on indi- 
rect distance measurement techniques, such as photometric 
parallax or main-sequence (MS) fitting (e.g.. lJohnsonlll95"7l : 
lEggen &"Sanda ge 1959), to probe the spa tial and kinemat- 
ical (sub-)structur es in the Galaxy (e.g., lJuric et af] 120081 : 
llvezic et alJl2008ah . 

Purely empirical photometric parallax relations can be con- 
structed based on trigonometric parallaxes to nearby stars. 
However, the difficulty arises from a severely restricted sam- 
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pie of stars with good trigonometric parallaxes over a wide 
range of stellar mass and metallicity. In particular, colors and 
magnitudes of MS stars are sensitive to chemical content. Un- 
fortunately, within the limited volume sampled by the Hippar- 
cos mission, reliable trigonometric parallaxes are available for 
only a handful of metal-poor stars ([Fe/H] < -LO), posing a 
challenge to Galactic halo studies relying on this information. 

A widely used technique for obtaining distances to indi- 
vidual stars and star clusters is to rely on theoretical predic- 
tions of stellar colors and magnitudes, tested with available 
observational constraints. Stellar evolutionary models pre- 
dict luminosities and effective temperatures (Teff) as a func- 
tion of time and chemical content (usually parameterized with 
[Fe/H], [a/Fe], and helium abundance Y). Stellar atmosphere 
models are then used to transform theoretical quantities (lumi- 
nosity and Teff) to observables such as magnitudes and colors, 
expressed in terms of bolometric corrections and color- Teff 
relations. Stellar evolutionary models can be tested against 
the Sun, nearby stars with accurate parallaxes, and other 
stars, such as eclipsing binaries that have accurate masses and 
radii. Furthermore, multicolor photometry in nearby clus- 
ters and field stars can be used to test bolometric correc- 
tions and color- Teff relations (e.g ., Vande nBerg & Clemli2003l : 
iPinsonneault et al.ll2004l:[An et aL,2007allbh . 

iPinsonneault et al] (I2003L l2004l) assessed the accuracy of 
theoretical stellar isochrones and reduced systematic errors 
in the model computation, particularly those arising from the 
transformation of theoretical to observational quantities. They 
demonstrated that stellar models from the Yale Rotating Evo- 
lutionary Code (YREC; ISills et alJlIOOOt) are in good agree- 
ment with observed mass-luminosity-Tetf relations for indi- 
vidual members (4500 < 7'eff(K) < 6000) in the Hyades open 
cluster with accurate Hipparcos parallaxes jde Bruiine et alj 
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1200 lb . However, they found that none of the widely-used 
color- Teff relations in the Johnson-Cousins-2MASS broad- 
band system (e.g., l Alonso et al.l I1995L 119961 : iLeieune et alj 
[1997, 1998) could reproduce the observed shapes of the MS 
on color-magnitude diagrams (CMDs) for the Hyades open 
cluster Since this implies problems with the adopted color- 
Jeff relations, they introduced e mpirical correction s to the 
color- Teff relations in the models. lAn et al.l (l2007bh showed 
that these corrections improve the shape match to MSs of 
other open clusters, as well as the internal consistency of dis- 
tances from several color indices. 

In this paper we perform a test of theo retical isochrones 
in the Sloan Digital Sky Survey (SDSS; 'York et al.' '2000^ 
[Stoughton et al. 2002; Abazaiian et al. 2003, 2004, 2005; 
lAdelman-McCarthy et al. 2006. 12007L 120081) ugriz photomet- 
ric system. Among previous and ongoing optical surveys, 
SDSS is the largest and most homogeneous database of stellar 
brightnesses currently available. SDSS measures the bright- 
nesses of stars using a dedicated 2.5-m telescope (iGunn et aH 
I2OO6) in five broadband filters u, g, r, /, and z, with av- 
erage wavelengths of 3551 A, 4686 A, 6 165 A, 7481 A, and 
8931 A, respectively (iFukugitaet all Il996t IStoughton eTaP 
I2OO2'). The imaging is carried out on moonless nights of 
good seeing (better than 1.6") under photometric conditions 
dHogg et alj|200l]) . A portion of the sky (along g reat circles) 
is irn aged in each run by 6 columns of CCDs (iGunn et alJ 
Il998l) . Astrometric positions i n SDSS are accu rate to < 0.1" 
for sources with r < 20.5 mag ( iPier et alj2003 ^. However, the 
SDSS filters were not originally designed for stellar observa- 
tions, owing to the primarily extragalactic mission of SDSS I. 
Therefore, it is important to understand the properties of stars 
in this system to exploit the full capabilities of the SDSS data 
set, especially for Galactic structure studies that require accu- 
rate distance measurements. 

Galactic globular and open clusters provide an ideal op- 
portunity to achieve this goal, because the same distance can 
be assumed for cluster members with a wide range of stellar 
masses. However, the stan dard SDSS photometric pipeline 
(Photo: ILupton et ani2002l) failed to provide photometry for 
the most crowded regions of high-density clusters. Photo was 
originally designed to handle high Galactic latitude fields with 
relatively low densities of field stars, and its photometry alone 
does not provide sufficiently well-de fined cluster se quences. 

In the previous paper of this series (lAn et al.l2008l hereafter 
Paper I), we employed the DAOPHOT/ALLFRAME (Stetson| 
[T987i, .199 4) suite of programs to derive photometry for 17 
globular clusters and 3 open clusters observed with SDSS. 
Our DAOPHOT photometry provides well-defined cluster se- 
quences from the lower MS to the red giant branch (RGB). 
The DAOPHOT photometry is on the native SDSS 2.5-meter 
system, so cluster fiducial sequences can be directly applied to 
other stars observed in SDSS, without relying upon any trans- 
formations. This will be of particular value for the study of 
the space motions of field stars in SDSS, where accurate dis- 
tances are required in order to make full use of the available 
proper motions. 

In light of this observational improvement, we assess the 
accuracy of theoretical isochrones that have been specifi- 
cally constr ucted in the SDSS 2.5-meter photometric system. 
IClem I (I2OO 6) compared isochrones to the fiducial sequences 
for three globular clusters (M13, M71, and M92) and two 
open cluster (M67 and N GC 6791), in the primed (u's'r'i'z! ) 
filter system (see also iRider et al.1 12004 iFornal et al.li2007h . 



However, the u'g'r'i'z' system defined by the ISmith et al] 
(l2002h sample of standard stars is different from the natu- 
ral ugriz sy stem of the SDSS 2.5-m survey telescope (see 
lAbazaiiane t al. 200 3i). Therefore, o ne must rely on transfor- 
mation equations (ITucker et al. 2006i), which limit the preci- 
sion of the derived ugriz magnitudes. 

We use several globular and open clusters that have been 
the subject of extensive studies in the literature. Distances 
to most of these clusters have been estimated bas ed on pho- 
tomet ry in the Johnson-Cousins system (e.g., ReidI 1 19971 
1998; Gratton et al. 1997; Carretta et al. 200Q, hereafter COO). 
Metal abundances for these clusters have also been obtained 
from high-resolution spectrosco py. In parti cular, a recent 
spectroscopic abundance study of iKraft & Iva ns (2003, here- 
after KI03) provides an opportunity to check for possible 
systematic errors in cluster distance and age estimates from 
the cluster abundance sca le, by comparison with those from 
ICarretta & GrattonI d 19971 hereafter CG97). 

As a base case, we test s tellar evolutionary models em- 
ployed in the previous work ( Pinsonneault et al]l2003L 120041 : 
lAn et al. 2007a.b) with updated input physics. We further 
insp ect theoretical isochrones from two independent stud- 
ies (iGirardi et al J 120041: lOotter et al.ll200 8>) to check theoret- 
ical uncertainties in the stellar interiors and atmosphere mod- 
els. Theoretical color-Tgff and bolometric corrections are 
als o compared with those from different atmosph ere mod- 
els (iHauschildt et"anil999l:ICasteni & Kuruc3i2()()3h to check 
theoretical uncertainties in the model atmospheres only. 

Our goal in this study is to determine whether these mod- 
els predict observed colors and magnitudes over a wide range 
of metallicity. In § 2 we summarize the photometric accuracy 
of the DAOPHOT/ALLFRAME cluster photometry in Paper I 
and the cluster properties (distance, reddening, and metallic- 
ity) adopted in this study. In § 3 we present a new set of 
stellar isochrones in ugriz and derive extinction coefficients. 
We al so compare our models with those from Girardi et al] 
(12004 and iDotter et all (12008). In § 4 we test the accuracy 
of these models using cluster fiducial sequences for MS stars, 
and present empirical corrections on color- Jeff relations for 
solar and super-solar metallicity clusters. In § 5 we derive 
MS-fitting distances and turn-off ages for our sample globular 
clusters, and demonstrate the ability to derive these quantities 
from the SDSS photometric database, using our set of stellar 
isochrones. 

2. THE CLUSTER DATA 
2.1. The Cluster Sample 

For the calibration of photometric parallax relations in any 
photometric filter system, it is important to select benchmark 
stellar systems with well-established stellar properties. In par- 
ticular, cluster distance, reddening, and metallicity estimates 
should be obtained reliably, independently from the photo- 
metric filter system in question. Although a few nearby open 
clusters have accurate geometric distance measurements (e.g., 
the Hyades and the Pleiades) in addition to precise reddening 
and metal abundance estimates, most of the stars in these clus- 
ters are too bright in SDSS because of the saturation limit in 
the survey (r « 14 mag). 

In this study, we selected six globular clusters (M3, M5, 
Ml 3, Ml 5, M71, and M92) and two open clusters (M67 
and NGC 6791) to test and calibrate theoretical stellar mod- 
els. These clusters were included in the sample of the 
DAOPHOT/ALLFRAME photometric data reduction in Pa- 
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TABLE 1 

Sample Globular Clusters and [Fe/H] 
Estimates in the Literature 



Name 


NGC 


Harris 


Zirm & West 


CG97^ 


KI03 


M15 


7078 


-2.26 


-2.15 


-2.15 


-2.42 


M92 


6341 


-2.28 


-2.24 


-2.15 


-2.38 


M13 


6205 


-L54 


-1.65 


-1.41 


-1.60 


M3 


5272 


-L57 


-1.66 


-1.34 


-1.50 


M5 


5904 


-L27 


-1.40 


-1.10 


-1.26 


M71 


6838 


-0.73 


-0.58 


-0.70 


-0.81 



^ Values are those shown in Table 1 of COO. 

per I. For our sample open clusters, there exist reliable dis- 
tances and metallicity estimates from high-resolution spec- 
troscopy, as discussed in more detail below. 

The MS of many globular clusters in our sample are re- 
solved --2-3 mag below the cluster MS turn-off (MSTO). 
They have spectroscopic metallicity estimates from CG97 and 
KI03, which are necessary to infer the absolute magnitudes of 
stars more accurately. Furthermore, distances to these clus- 
ters have been derived in the literature based on Hipparcos 
parallaxes to nearby subdwarfs, either in the Johnson-Cousins 
or Stromgren photometric systems. We therefore view these 
clusters as ideal tests of the color transformations for SDSS 
filters. We first review the DAOPHOT photometry in Paper I, 
and discuss cluster properties adopted in the model compari- 
son (§11. 

2.2. Cluster Photometry and Fiducial Sequences 

Cluster photometry was obtained in Paper I from the SDSS 
ugriz imaging data, employing the DAOPHOT/ ALLFRAME 
suite of programs". Due to the limited contamination from 
background stars, fiducial sequences of our sample clusters 
could be accurately derived on CMDs with several color in- 
dices (u- g, g- r, g-i, and g-z) and r as a luminosity index: 
hereafter, {u-g,r), {g-r,r), (g-i,r), and {g-z,r), respec- 
tively. 

The zero-point error of the DAOPHOT photometry is ^ 
l%-2%, estimated from repeated flux measurements of stars 
in overlapping SDSS strips/runs. The DAOPHOT magnitudes 
were tied to the Photo magnitude system using a set of clus- 
ter flanking fields, which are far away from the dense cluster 
cores. Error distributions at the bright ends indicate photomet- 
ric r andom error s of ~ 1% in griz and ~ 2% in the u band: see 
also'Ivezic elal] (|2004|) . 

In Paper I we found that the zero points of the DAOPHOT 
photometry for M71 are uncertain, because of inaccurate 
Photo magnitudes in the cluster's flanking fields. For this 
reason, we adopt ed a fiducial sequence in u'g'r'i'z' from 
IClem et al] (1200 8'). We included this cluster in the following 
analysis because M71 has an intermediate metal abundance 
([Fe/H] « -0.8), bridging the gap between globular and open 
clusters. We converted the Clem et al. fiducial sequences in 
the u'g'r'i'z' system onto the S PSS 2.5-m ugriz sys tem, using 
the transformation equations in lTucker et al.l (l2006h . 



2.3. Cluster Properties: Globular Clusters 

Table [U lists [Fe/H] estimates of clusters in the literature. 
In the model comparison, we consider metallicity scales for 
globular clusters from CG97 and KI03 (see also Kraft & 
Ivans 2004). In each study a consistent technique was em- 
ployed to derive metallicities for giants in all of our sample 



clusters. However, the CG97 values are based on Fe I lines, 
while the KI03 [Fe/H] values are based on Fe II lines from 
high-resolution spectra, which were argued to be less affected 
by non-local thermodynamic equilibrium (non-LTE) effects. 
KI03 found that their [Fe/H] values are about 0.2 dex lower 
than those of CG97 and that they agree well with those of 
IZinn & Wesll d 1 9841) for intermediate metallicities. 

Table |2] lists distances for the globular clusters in our sam- 
pl e. The secon d column shows distances from the catalog 
of lHarrisI d 19961) February 2003 revision, which were primar- 
ily based on the mean V magnitude of the horizontal branch 
(HB). The next five columns list distances derived from the 
Hipparcos-hased subdwarf fitting in the literature (Reid 1997, 
1998; COO; Grundahl et al. 2002; KI03). These studies used 
a subset of the Hipparcos subdwarfs to fit the MS of a clus- 
ter on CMDs, in either the Johnson-Cousins or th e Stromgren 
photometric system. The quoted precision from iReidI (1 19971 
Il998i) is shown in the table. We computed true distance mod- 
uli of clusters in COO, taking their adopted ratio of total to 
selective extinction Ry = Av/E(B-V) = 3.1 and E(B-V). 
The errors are the sum (in quadrature) of their fitting and to- 
tal systematic errors ( 0.12 mag). The s ixth c olumn shows the 
distance to M71 from lGrundahl et aP (l2002h . who employed 
Stromgren uvby ph otometry for b oth subdwarfs and the clus- 
ter in MS fitting. iGrundahl et a l.' found an apparent distance 
modulus (m-M)y = 13.71 ± 0.04(ran) ± 0.07(sys) for their 
adopted reddening of E{B-V) = 0.275, which is translated 
into (m-M)o = 12.86 ±0.08. KI03 used the Hipparcos sub- 
dwarfs to derive distances to five globular clusters from their 
application of the MS fitting technique. The last two columns 
in Table 12] list distances from MS fitting in this study, which 
are discussed in § 15.11 

All of the above studies were based on a similar set of Hip- 
parcos subdwarfs. However, there exist about a 10% dif- 
ference in their distance estimates. This is because several 
systematic errors were involved in the subdwarf-fitting tech- 
nique, such as the metallicity scale for subdwarfs and clus- 
ters, foreground reddening, effects of binaries, and the pho- 
tometry used. Since no study has provided homogeneous 
distances to all of our sample clusters, we assumed the fol- 
lowing subdwarf-fitting distances in the model comparison 
(§11: COO (M5 M13 , and M92), KI03 (M3 and M15), and 
IGrundahl et aP dlOOl M7 1). 

We adopted the COO distance estimates for M5, M13, and 
M92, since they employed the same metallicity scale for both 
subdwarfs and cluster stars in the subdwarf-fitting technique. 
Nevertheles s , distances to M5 and Ml 3 in COO agree with 
iReidI (ll997h . lReidl(ll998h . and KI03, within their specified er- 
rors, although these studies adopted slightly different metal- 
licity scales for both subdwarfs and clusters. Notable differ- 
ences can be found for M92, where the Reid (1997) value 
is about 0.3 mag larger than the COO estimate, and about 
0.2 mag larger than the KI03 value. As discussed in COO, 
a part of this discrepancy can be explained by the inhomoge- 
neous metallicity scale used in lReidI (Il997h at the low metal- 
licity end. For this reason , we also adop t ed the distance to 
M15 from KI03. Although GrundahleLaij dlOOl utilized the 
Stromgren uvby photometry to estimate the cluster metallicity, 
rather than high-resolution spectroscopy, they provided an in- 
ternally more consistent metallicity scale for both subdwarfs 
and cluster stars than Reid (1998); see Figures 10 and 11 in 
IGrun dahl et aPdlOOl) . 

Table [3] summarizes reddening estimates for our sample 



" Available at jhttp : //www, sdss . org/dr6/products/value_addedgk Mjlj>kgsCjWEeiSlifttfkep^ .fcflft columns list values 
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TABLE 2 

Distance Moduli of Globular Clusters 



Hipparcos-haied Subdwarf Fitting 



Tliis Study 



Cluster Harris" Reid(1997) Reid(1998) 



COO" 



Grundahl et al. KI03 [Fe/H],- 



[Fe/H] 



KI03 



M15 

M92 

M13 

M3 

M5 

M71 



15.06 
14.58 
14.42 
15.09 
14.37 
13.02 



15.38±0.10 
14.93 ±0.10 
14.48±0.10 

14.45 ±0.10 



14.52±0.15 
13.19±0.15 



14.64 ±0.14 
14.38±0.13 

14.46 ±0.1 3 



15.25 
14.75 
14.42 
15.02 
14.42 



12.86±0.08 



15.14±0.11 
14.63 ±0.03 
14.41 ±0.05 
15.00 ±0.07 
14.35 ±0.07 
12.96 ±0.08 



15.10±0.11 
14.66 ±0.04 
14.34 ±0.05 
14.95 ±0.07 
14.26 ±0.06 
12.86 ±0.08 



" Distance moduli assuming Ay /E(B -V) = 3.1 and reddening values from the catalog of lHarri^tl99(i) February 2003 r 
of systematic and random enors. 



Quadrature sum 



TABLE 3 
Previous £"(5 -V) Estimates 
for Globular Clusters 



Cluster 


Hanis 


Schleeel et al." 


KI03 


M15 


0.10 


0.11 


0.10 


M92 


0.02 


0.02 


0.02 


M13 


0.02 


0.02 


0.02 


M3 


0.01 


0.01 


0.01 


M5 


0.03 


0.04 


0.03 


M71 


0.25 


0.32 


0.32 



" Quoted precision (16%). 



taken from the Harris compilation, the 'Schl egel et alJ (Il998h 
dust maps, and E(B-V) estimates from KI03, who compared 
colors derived from high-resolution spectroscopic determina- 
tions of Tgff with the observed colors of the same stars. There 
is no appreciable offset in the mean between the three sets of 
estimates. 

Since distances from MS fitting are sensitive to reddening, 
the same reddening values as in the individual subdwarf fitting 
studies should be used in the model comparison. However, we 
adopted E(B-V) estimates from KI03 for all of our sample 
globular clusters, because they provided homogeneous stellar 
E{B-V) estimates from a consistent technique. Neverthe- 
less, the E{B-V) values adopted in COO are only 0.005 mag 
larger than those in KI03 for M5 and M92, with the same 
E{B-V) for M13. In addition, with the exception of M71, 
the difference between th e KI03 values and those of the dust 
map values from Schleg el et al.l ([T998) is 0.004 ±0.005 mag 
in the sense of lSchlegeletal.1 (Il998h minus KI03. In the 
case of M7I, we ado pted the same E{B -V) = 0.28 as in 
iGrundahl etal] (l2002h . We assumed a 20% eiTor in E(B-V) 
in the model comparison (§ |4|i for all clusters in our sam- 
ple, which encompasses the differences between the studies 
shown in Table [3] 

2.4. Cluster Properties: Open Clusters 

Table |4] lists metallicity, reddening, and distance estimates 
for our sample open clusters. For M67 we took the aver- 
age reddening and metallicity estimates from high-resolution 
spectroscopy in the literature (l An et alJl2007bl and references 
therein). Although M67 is closer than any other cluster in this 
study, it is still too far away for Hipparcos to provide a direct 
distance measurement. We adopted a cluster distance esti- 
mated from an empirically calibrat ed set of isochrq nes in the 
Johnson-Cousins-2MASS system (I An et al.ll2007b l). In that 
paper, we noted that the application of these isochrones re- 
sulted in consistent estimates for metallicity and reddening 
with those in the literature. 

For NGC 6791 there exist four recent [Fe/H] measure- 



ments from high-resolution spectroscopic stud ies: [Fe/H] = 
+0.47 ± 0.04 (random) ± 0.0 8 (systematic) jGratton et alj 
l2006l) . +0.39 ±0.01 (random) dCmraro et al.^'2006'), +0.35 ± 
0.02 (r andom) ± 0.1 (sys tematic) (Origlia et al. 2006), and 
+0.40 (Jen sen et a"Lll2006l) . In the model comparison (§ g]), 
we adopted the average of these measurements ([Fe/H]) = 
+0.40 ± 0.03 with the same weight given to each study. Us- 
ing an extended set of calibrated isochrones in the Johnson- 
Cousins-2MASS system, we derived a photometric metal- 
licity, [Fe/H] = +0.43 ±0.07 (M. H. Pinsonneault & D. An 
2009, in preparation), that is consistent with the above mea- 
surements, assuming solar abundance mixtures (see Gratton 
et al. 2006 and Origlia et al. 2006 for subsolai- [C/Fe]). How- 
ever, we found E{B-V) = 0.098 ±0.014 from the calibrated 
isochrones, while reddening estimates in the literature vary 
between 0.1 < E(B-V) < 0.2. We adopted our E(B-V) 
in the model comparison. Both our photometric metallic- 
ity and reddening values yield a cluster distance modulus, 
(m-M)() = 13.02 ±0.05, assuming solar abundance mixtures. 
If we use our adopted metallicity [Fe/H] = +0.40 in this study, 
the distance modulus decreases by « 0.03 mag. Given the size 
of the Icr error in distance, however, we neglected this dif- 
ference, and adopted (m-M)o = 13.02 ±0.05 for the model 
comparison. 

3. MODELS 

In this section we generate a set of theoretical isochrones 
in ugriz, which serves as a base case for the comparison with 
cluster fiducial sequences (§ ID) and in the estimation of dis- 
tances and ages of globular clusters (§|5]l. We further use other 
theoretical isochrones and color- Teff relations in the literature 
to inspect theoretical uncertainties in the model computations. 

3.1. Theoretical Isochrones in ugriz 

Stellar e v olutio n ary tracks were generated using YREC 
dSills et all l200tf , 'Dela haye & Pinsonneault! |2006') for a 
wide range of compositions and ages. Similar models 
were employed in our previous studies in the Johnson- 
Cous i ns-2MASS filt er system ( Pinsonneault et alj I2003L 
l2004tlAn et aPIIOOTallbl) . but the new models include updated 
input physics (see below). Details of the interiors models will 
be presented in a paper in preparation. 

We recently adopt ed atomic opa city data from the Opac- 
ity Project (OP, Badnell et al.ll2005)) for our models. Prev ious 
models used the OPAL opacities (llglesias &Rogersll996l) . In 
our view the OP data have benefited from improvements in 
the equation of state, especially for iron, and the availabil- 
ity of the monochromatic opacities permits tables to be read- 
ily constructed for a wid er range of mixtures than for OPAL. 
See lBadnelleTan ( l2005h for a comparison of the underlying 
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TABLE 4 
Properties of Open Clusters 





Alternate 


High Res. 


E(B 


-V) 


(m-M)o 




NGC 


Name 


[Fe/H] 


Previous Work 


Schlegel et al." 


Previous Work This Study 


References 


2682 


M67 


+0.00 ±0.01 


0.041 ±0.004 


0.032 


9.61 ±0.03 9.59 ±0.03 


1 


6791 




+0.40 ±0.03 


0.098 ±0.014 


0.155 


13.02 ±0.05 12.95 ±0.05 


2 



References. — References for high resolution spectroscopic [Fe/H] values in the literature, and estimates for stellar E{B - V) and 
(m— M)o: (1) An et al. 2007b, and references therein; (2) Pinsonneault & An (2008, in preparation), and references therein. 
^ Quoted precision (16%). 



opacities and a discussion of the differences in the underly- 
ing physica l model. For T < 10"* K we used the molecular 
opacities of [Alexander & Fergusonl(ll994h . We recognize that 
there is an inconsistency between the mixture of heavy ele- 
ments in these tables and the Grevesse & Sauval ( 1998) mix 
adopted in this paper, but the mixtures are quite similar, and 
the impact on the opacities is typically at the 0.01-0.03 dex 
level for solar metallicity. We used the OPAL 2001 equation 
of stat e ("Rogers & Navfo novl 12002 ) for T > 10^ K and the 
ISaumo n et al. ( 1995) equation of state for T < 10^^ K. In the 
transition region between these two temperatures, both for- 
mulations are weighted with a ramp function and averaged. 
The chemical compositi on of each shell wa s updated using 
nuclear reaction rates in Bahcall et al.l (1200 ll) . The models do 
not include gravitational settling. A mixing-length of a = 1 .78 
and the solar helium abundance Fq = 0.269 were calibrated 
by matching the solar radius (Rq = 6.9598 x 10** m) and lu- 
minosity (Lq = 3.8418 X 10^^ ergs s~') at t he age of the Sun 
(4.57 Gyr); see iDelahave & Pinsonneaulj (|2006|) . We used 
the ATLAS9 atmosphere models to set the surface boundary 
condition, defined as the pressure at r = 2/3. 

At [a/Fe1 = 0, we used the solar mixture of 

iGrevesse & Sauval (Il998i) for the initial chemical mixture, 
and generated models at [Fe/H] = -3.0, from [Fe/H] = -2.0 
to [Fe/H] = -0.5 with intervals of A[Fe/H] = 0.25, 
from [Fe/H] = -0.5 to [Fe/H] = 0.5 with intervals of 
A[Fe/H] = 0.1, and [Fe/H] = 0.75. We assumed the 
helium enrichment parameter, AF/AZ = 1.2, along with 
the primordial helium abundance, Yp = 0.245 ± 0.002 (see 
iPinsonneault et all 120031) '^. We also computed a-element 
enhanced models at [Q:/Fe]= +0.3 and +0.6 for each [Fe/H], 
by increasing the ratio of a elements (O, Ne, Mg, Si, S, Ar, 
Ca, an d Ti) by the same amount in the IGrevesse & Sauval 
( Il998l) mixture using OP tables computed with the the same 
mixture. We followed the evolution of each model from 
pre-MS until hydrogen in the core is exhausted. A discussion 
of the post-MS evolutionary tracks is outside the scope of this 
paper 

Theoretical isochrones in the luminosity- Teff plane were 
then transformed onto the observed C MP plane, using a set of 
synthetic specti-a. We used MARCS (iGustafsson et al.ll2008l) 
to construct stellar atmosphere models over a wide range of 
Teff (from 4000 K to 8000 K with intervals of 250 K), logg 
(from 3.5 dex to 5.5 dex with intervals of 0.5 dex), and chem- 
ical abundances (see below). These models are based on 
1-D plane-parallel geometry, hydrostatic equilibrium, mix- 
ing length theory (MLT) convection, and LTE. Microturbu- 
lence (^) was chosen to be 2.0 km s~' . We assumed the same 
[a/Fe] abundance scales (IGrevesse & Sauvailll99 8') as those 
in the interior model computations. Based on these mod- 

We assumed that the helium abundance is a function of the heavy- 
element content Z as given by Y = Yp + ( AY / AZ)Z. 



els, we constructed synthetic spectra at [Fe/H] ([a/Fe]) = 
-3.0 (0.4), -2.0 (0.0,0.4), -1.5 (0.0,0.4), -1.0 (0.0,0.4), 
-0.5 (0.0,0.2), -0.3 (0.0,0.1), -0.2 (0.0,0.1), -0.1 (0.0), 
0.0 (0.0), +0.1 (0.0), +0.2 (0.0), +0.4 (0.0). 

Stellar magnitudes were derived by convolving synthetic 
spectra with the SDSS ugriz filter response curves'^, which 
include extinction through an air mass of 1 .3 at Apache Point 
Observator y. We integrated flu x with weights given by pho- 
ton counts (Gir ardi et al .1120021) . Magnitudes were then put 
onto the AB m agnit ude system using a flat 363 1 Jy spectrum 
(Oke & Gunnll983l) . Bolometric corrections in the r band and 
synthetic colors are presented in Table |5] where we adopted 
the bolometric ma gnitude of the Sun, Mboi © = 4.74 mag 
dBessellet al.lll998l) . 

In the following analysis, we adopted an a-element 
enhancement scheme motivated by the observational na- 
ture of these elements among field and cluster stars from 
high-resolu t ion s pectr oscopic stud i es, a s summarized by 
IVenn et al.1 (|2004 and iKirbvet al.1 (l2008l) : [a/Fe]= +0.4 at 
[Fe/H]=-3.0, [a/Fe]=+0.3 at [Fe/H]= -2.0, -1.5, -1.0, 
[a/Fe]= +0.2 at [Fe/H]= -0.5, and [a/Fe]= +0.0 at [Fe/H]= 
-0.3, -0.2, -0.1, +0.0, +0.1, +0.2, +0.4. We used a linear in- 
terpolation in this metallicity grid to obtain isochrones at an 
intermediate [Fe/H] value. As an alternate case, we also tested 
models with [afFe]=+OA at [Fe/H] = -2.0,-1.5,-1.0, given 
a A[Q:/Fe] ^ ±0.1 dispersion in the observed distribution of 
[a/Fe]. 

Isochrones constructed in this way are on a perfect AB 
magnitude system, in which magnitudes can be translated di- 
rectly into physical flux units. However, it is known that the 
SDSS photometry slightly deviates from a true AB system 
(fAbazajian et al. 2004). To compare our models with fidu- 
cial sequences, we adjusted model magnit udes using AB cor- 
rections given by lEisenstein et al.l (l2006l) : mab = m - 0.040, 
'AB = / + 0.015, and zab = z + 0.030, with no corrections in 
g and r (see also Holberg & Bergeron 2006). 

T he SDSS magnitud es are on the asinh magnitude sys- 
tem (iLupton et al.ll 19991) . which is essentially the same as the 
standard astronomical magnitude at high signal-to-noise ra- 
tios.''* Near the 95% detection repeatability limit for point 
sources in SDSS (u = 22.0 mag, g = 22.2 mag, r = 22.2 mag, 
/ = 21.3 mag, and z = 20.5 mag), the difference between the 
asinh and Pogson magnitudes are < 0.01 mag, and it rapidly 
becomes zero as a source becomes brighter. Since our fiducial 
sequences were derived using relatively bright stars (Paper I), 
we compare theoretical isochrones on the Pogson system di- 
rectly with the DAOPHOT cluster photometry on the asinh 
system, without relying upon any transformation. 

The isochrones constructed in this way are available at 
|http : / /spider . ipac . caltech . edu/ ~deokkeun/sdss . htni 

^'^ see also |httpT77wwwT^dssT^rg7dr67aTgorTthm^ ■ html| 



6 



An et al. 



TABLE 5 

BoLOMETRic Corrections and Synthetic Colors in ugriz. from the 
MARCS Model Atmospheres 



[Fe/H] 


[a/Fe] 


Tcff 




BC(r) 


u-g 


g-r 


g-' 


g-z 


(dex) 


(dex) 


(K) 


(dex) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


-3.0 


0.4 


4000 


3.0 


-0.398 


2.037 


1.115 


1.584 


1.867 


-3.0 


0.4 


4000 


3.5 


-0.374 


1.988 


1.111 


1.579 


1.859 


-3.0 


0.4 


4000 


4.0 


-0.354 


1.967 


1.116 


1.583 


1.860 


-3.0 


0.4 


4000 


4.5 


-0.333 


1.973 


1.127 


1.592 


1.866 


-3.0 


0.4 


4000 


5.0 


-0.306 


1.999 


1.146 


1.609 


1.878 


-3.0 


0.4 


4000 


5.5 


-0.275 


2.057 


1.169 


1.629 


1.894 



Note. — The AB con'ectioiis for the SDSS magnitudes and empirical color conections are 
not applied. Table 5 is available in its entirety in the electronic edition of the Asrwphysical 
Journal. A portion is shown here for guidance regarding its form and content. 



In the following analysis, we refer to our models as 
YREC+MARCS to dist ingu ish from the other theoretical 
isochrones described in § 13.31 



3.2. Extinction Coefficients 

We derived extinction coefficients in the ugriz filters us- 
ing MARCS synthetic spec tra. We first r eddened the syn- 
thetic spectra assuming the ICardell i et alJ (fi989) extinction 
law, which is given in polynomial equations for the optical 
wavelength range as a function of the wavelength and Ry- We 
then convolved the reddened spectra with the SDSS filter re- 
sponses, over a wide range of stellar parameters and Ry- In 
addition to the ugriz filter response functio ns, we also used 
the normalized V filter response fromlBessell (1990) to com- 
pute the absorption ratios of SDSS filters with respect to the 
Johnson V band, Ajy/Ay, where A represents one of the SDSS 
filters. 

For solar metallicity models with Teff = 5750 K and logg = 
4.5, we foundA„/Ai/= 1.567 ±0.044, Ag/Ay = 1.196±0.013, 
Ar/Av = 0.874 ± 0.005, A/ /Ay = 0.672 ±0.011, and A, /Ay = 
0.488 ±0.014, assuming Ry = 3.1. The errors represent the 
case of assuming ARy = ±0.2. Although there exists a mild 
variation as a function of colors or Teff of models, differences 
are within « Icr of the above estimates. In this paper we ne- 
glected these color-dependent variations in extinction ratios, 
since most of the sample clusters have interstellar reddening 
less thangCfl- V) ~q .l. 

Girardi et alj ( |2004|) also estimat ed extinction coeffic ients 
based on the same reddening law (ICardelli et al.lll989 ). but 
using synthetic ATLAS9 spectra. Our values are in good 
agreement with their values. At 4000 < 7;ff(K) < 8000 with 
logg = 4.5, the differences are only AA^/Ay ~ 0.01. 

We additionally derived extinctio n coefficients in u 'g'r'i'z' 
for the fiducial sequence of M71 (IClem et al.l l2008-). since 
transformations fr om ugriz to u'g'r 'i'z' were derived for 
unreddened colors (iTucker et al.l2006l) . We derived extinction 
coefficients for USNO u'g'r'i'z' filters'^ using solar metallic- 
ity models with Teff = 5750 K and logg = 4.5 at Ry = 3.1: 
A„,/Ay = 1.569, Ag,/Ay = 1.177, A^'/Ay = 0.866, A,/ /Ay = 
0.656, and A,VAv= 0.488. 

For a given Ry and ^(B-y) of a cluster, we computed ex- 
tinction in the r band and color-excess values in ugriz from 



and 



A,= {^jEiB-V)Ry, (1) 



E{A,-A2)={^-^]E{B-V)Ry, (2) 

Ay Ay 



See |http : //www-star . f nal ■ gov/ugriz/Filters/response" 



where Ai and A2 represent each of the SDSS filter passbands. 

3.3. Other Models: Isochrones 

Isochrones in ugriz are also availa ble from two other inde- 
pendent studies. i Girardi et al.l ([2004) provided the first exten- 
sive set of isochrones in ugriz filters, b ased on evolutionar y 
models used in their earlier studies (see iGirardi et al.ll2002l) . 
To derive bolometric corrections and color-Teff relations, 
they have mainly used ATLAS9 non-o vershooting models 
('Castelli et al.1 119971: iBessell et al] |1998): hereafter we refer 
to these isochrones as PadovaH-ATLAS9. On the other hand, 
iDotter et aP (l2008l) recently presented the Dartmouth Stel- 
lar Evolution Progr am (DSEP) isochrones, where PHOENIX 
model atmospheres (iHauschildt et al.l 19991) were employed to 
convert luminosity and Teff into magnitudes and colors; here- 
after we refer to these isochrones as DSEPh-PHOENIX. In this 
section we compare different evolutionary models (L - Teff) 
and atmosphere models (bolometric corrections and color- Teff 
relations) separately. 

Figure [T] shows the comparison between Padova (dotted 
lines) and YREC (solid) isochrones on the log(L/L0) - Teff 
plane. Padova isochrones are shown at Z = 0.0001 ([Fe/H] = 
-2.32), Z = 0.0004 ([Fe/H] = -1.71), Z = 0.001 ([Fe/H] = 
-1.31), Z = 0.004 ([Fe/H] = -0.71), Z = 0.019 ([Fe/H] = 
+0.00), and Z = 0.040 ([Fe/H] = +0.37), taking their adopted 
solar composition (Y = 0.273 and Z = 0.019) to compute 
[Fe/H] values. All of these models are ba sed on the solar 
abundance ratio, since IGirardi et al.l (|2004|) do not provide 
models with a-element enhancement. 

As shown in Figure [T] YREC and Padova isochrones are 
generally in good agreement with each other, but they exhibit 
a large difference at the lower MS, in the sense that the Padova 
models are hotter than the YREC models for a given lumi- 
nosity. This difference is likely due to the use of different 
equations of state at low temperatures. Padova models use 
those in Mihalas et al. (1990), while YREC models adopt the 
ISaumon et al.l (119951) equation of state. The morphologies of 
the MSTOs differ for the solar metallicity, 3.5 Gyr models, 
bec ause Padova models include convective core overshooting 
(see'Girardi^eLalJ|2000|). 

Figure |2] shows a similar comparison in the log(L/L0)- 
Teff plane between DSEP (dotted line) and YREC (solid 
line), both of which include a-element enhancement. DSEP 
isochrones at intermediate [a/Fe] values were obtained by in- 
terpolating models over their customized grid at [a/Fe]= 0.0, 
0.2, and 0.4. In the left-hand panel, these isochrones are 
shown at [Fe/H] =-2.42, -1.60, -1.26, and -0.81, with 
our adopted a-element enhancement scheme ([a/Fe]« +0.3 
lraiit Fe/H]< -1; see § 13. 11 1. In the right-hand panel, isochrones 
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Fig. 1. — Comparison between Padova (dotted line) and YREC (solid line) evolutionary models. Left: Isochrones are shown at [Fe/H] = —2.32, —1.71, —1.31, 
and— 0.71, with [a/Fe] = 0.0. Ages of 12.6 Gyr are assumed. Right: Isochrones are shown at [Fe/H] = 0.0 (3.5 Gyr) and +0.37 (10.0 Gyr) with [a/Fe] = 0.0. 



are shown at [Fe/H]= +0.0 and +0.4 with no a-element en- 
hancement. 

As shown in Figure |2] the two models are generally in 
good agreement with each other on the log(L/Lo) vs. T^ff 
plane. In fact, DSE P has been develop ed based on a simi- 
lar code to YREC (Dema rque et al.ll2007 ). and many input in- 
gredients are the same as in YREC. However, they have been 
subsequently modified by the two independent groups. For 
example, DSEP includes the PHOENIX model atmosphere 
to set the surface boundary condition, while YREC adopts 
the ATLAS9 model. DSEP uses high temperature opacities 
from OPAL (Iglesias & Rogers 1996), while YREC utiHzes 
OP data (B adnell et al.l2005l) . As mentioned above, the differ- 
ent MSTO morphology of 4 Gyr, solar-metallicity isochrones 
is due to the inclusion/exclusion of the convective core over- 
shoot in DSEP/YREC. 



3.4. Other Models: Bolometric Corrections and Color-T^s 
Relations 

Figure [3] compares r-band bolometric corrections and syn- 
thetic colors from ATLAS9"' and MARCS, in the sense of the 
former minus the latter. We derived magnitudes and colors for 
ATLAS9 spectra, following the same steps as for MARCS. 
Comparisons are shown at \ogg = 4.0 with three different 



metallicities: [Fe/H]= -2.0 and [a/Fe]= 0.4 (dotted line), 
[Fe/H]= -1 .0 and [a/Fe]= 0.4 (dot-dashed line), [Fe/H]= +0.0 
and [Q!/Fe]= 0.0 (solid line). Although the difference between 
ATLAS9 and MARCS is becoming larger at a lower T^ff end 
in the comparison, the agreement between these two models 
is generally good (< 0.02 mag in bolometric corrections and 
colors) at 5000 < Teff (K) < 8000. The comparison in u-g 
suggests a large theoretical uncertainty in colors that involve 
the M bandpass. 

Similarly, Figure |4] shows comparisons in r-band bolomet- 
ric corrections and synthetic colors between MARCS and 
the LTE-based NextGen grid generated with the PHOENIX 
code.'^ Comparisons are shown at logg = 4.0 with three dif- 
ferent metallicities: [Fe/H] = -2.0 (dotted line), -1.0 (dot- 
dashed line), and +0.0 (solid line). However, a-element en- 
hanced models were not included in this comparison, because 
they were not available in the public PHOENIX database. 
The differences are in the sense of PHOENIX/NextGen minus 
MARCS. Unlike ATLAS9 model colors, PHOENIX/NextGen 
models have larger color differences from MARCS in all color 
indices, especially in u-g. 

4. MODEL COMPARISON 

In this section, we assess the accuracy of theoretical 
isochrones using cluster fiducial sequences for our sample 



^ |http : / /wwwuser . oat .ts.ast ro. it /castelli/ grids. html| '^ |http : / /phoenix . ens-lyon . f r/ Grids/NextGen/ SPECTRA/) 
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Fig. 2. — Comparison between DSEP {dotted line) and YREC {solid line) evolutionary models. 
and -0.81, with ct-element enhancement adopted in this study ([a/Fe]Si +0.3 at [Fe/H]< -1; see §[371 
shown at [Fe/H] = 0.00 (3.5 Gyr) and +0.40 (10.0 Gyr) with [a/Fe] = 0.0. 



Isochrones are shown at [Fe/H] = -2.42, -1.60, -1.26, 
Ages of 12.6 Gyr are assumed. Right: Isochrones are 



clusters. For our base case, we use YREC+MARCS mod- 
els to check whether stellar models satisfactorily reproduce 
observed colors and magnitudes of stars over a wide range 
of metal abundance. In this comparison, we take into ac- 
count various systematic errors, such as those from the photo- 
metric calibration, cluster distance, and the metallicity scale 
for globular clusters. We also test PadovaH-ATLAS9 and 
DSEPh-PHOENIX isochrones with cluster fiducial sequences. 

4.1. Systematic Errors in Color-Ten Relations 



Figures |5]-[8] show comparisons between YRECh-MARCS 
models and fiducial sequences in (g-r,r), (g-i,r), (g-z,r), 
and {u-g,r), respectively. In each panel, fiducial sequences 
are displayed as either filled or open circles, and the mod- 
els are shown as solid lines with [Fe/H] as indicated on top of 
each panel. The leftmost panel contains fiducial sequences for 
the two most metal-poor clusters (Ml 5 and M92). The sec- 
ond and third panels display those for intermediate-metallicity 
globular clusters (M3, M5, M13, and M71). The last panel 
displays fiducial sequences for two open clusters (M67 and 
NGC 6791). To place fiducial sequences on dereddened color 
vs. absolute magnitude diagrams, we adopted the cluster dis- 
tance and reddening values as discussed in § |2l The range of 
M, covered by a fiducial sequence is different from one cluster 
to the other, because SDSS images were taken in drift-scan or 
time-delay-and-integrate (TDI) mode, with a fixed effective 



exposure time of 54. 1 seconds per source. For the models, we 
assumed ages of 12.6 Gyr for all globular clusters, 3.5 Gyr 
for M67, and 10.0 Gyr for NGC 6791. To guide the eye, the 
approximate positions of the Sun'^ are shown on the CMDs 
as bulls-eyes. 

The model comparisons in Figures [SHS] are reproduced in 
Figuresl9HT2las a function of T^s in g-r, g-i, g-z, and u-g, 
respectively. We estimated color differences between mod- 
els and cluster fiducial sequences at a given M^, which was 
then converted into T^ff in the model. Open circles represent 
color differences for each of the fiducial points, as shown in 
Figures |5}{8] In this comparison we used models at the KI03 
[Fe/H] values for the globular clusters. Dotted lines indicate 
zero differences between model and the data for each cluster 

In Figures IgUTSl the difference between open and filled cir- 
cles effectively show the size of a systematic error in the 
model comparison from photometric zero-point errors. In Pa- 
per I we used Photo magnitudes in a set of flanking fields, 
which are far from the crowded cluster field, to set the zero 
point of DAOPHOT cluster photometry. This procedure was 
necessary to put DAOPHOT cluster photometry on the Photo 
magnitude scale. There were typically two SDSS imaging 
runs covering each cluster, and DAOPHOT photometry was 
independently calibrated on each run. However, about 1 %- 

See http : 77www . sdss ■ org/dr 6 /algorithms /sdssUBVRI Trans form, h' 
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Fig. 3. — Comparisons in r-band bolometric corrections and synthetic 
colors between ATLAS9 and MARCS atmosphere models, in the sense of 
ATLAS9 minus MARCS. Comparisons are shown at logg = 4.0, with three 
different metallicities: [Fe/H]= -2.0 and [a/Fe]= 0.4 (dotted line), [Fe/H]= 
-1.0 and [a/Fe]= 0.4 (dot-dashed line), [Fe/H]= 4-0.0 and [a/Fe]= 0.0 (solid 
line). 



m 



I 



I 



I 

<] 



N 
I 



-0.10 
0.00 
0.10 

0.00 
0.10 
0.20 
0.30 
-0.10 
0.00 
0.10 

-0.10 
0.00 
0.10 

-0.10 
0.00 
0.10 



1 — I — \ — \ — \ — \ — I — \ — I — I — I — I — I — I — I — \ — I — \ — \ — \ — \ — l- 
- NextGen vs. MARCS 




I I I I I I 



I I I I I I 



_L 



_L 



_L 



_L 



_L 



8000 7000 6000 5000 4000 
Terr (K) 

Fig. 4. — Comparisons in r-band bolometric corrections and synthetic col- 
ors between PHOENIX/NextGen and MARCS atmosphere models, in the 
sense of PHOENIX/NextGen minus MARCS. Comparisons are shown at 
logg = 4.0, with three different metallicities: [Fe/H] = -2.0 (dotted line), 
—1.0 (dot-dashed line), and 4-0.0 (.mlid line). All models assume the solar 
abundance ratio. 



2% systematic differences were found in DAOPHOT photom- 
etry for overlapping stars in these two runs. For this reason, 
we combined photometry from two runs by selecting a "ref- 
erence" run, which was simply chosen based on the sky cov- 
erage of a cluster or the run number. Open circles in Fig- 
ures I9HT2I and the data in Figures |5]-[8] show the case when 
we adopt these zero points for our sample clusters. On the 
other hand, filled circles show cases when using an alternate 
set of zero points for the same stars. In most instances, this 
is the zero point established by the other run. For M67 and 
NGC 6791, there were images taken in three SDSS runs, but 
we only display the case that shows the largest difference in 
magnitu des. In M71 , open c ircles show the model comparison 
with the IClem et al] (|2008|) fiducial sequence (after a photo- 
metric transformation to ugriz), where we adopted a 0.01 mag 
zero-point error in each passband (see below). 

The grey strip in Figures l9l-[T2l represents a ±lcr range in 
the color comparison from various systematic errors. We as- 
sumed the difference between open and filled circles as a 2a 
error from the photometric calibration. We adopted quoted er- 
rors in distance from individual studies (§|2]). Since KI03 did 
not provide the distance error, we assumed (T(,„_m) = 0.15 mag, 
which is a reasonable size for the errors of the Hipparcos- 
based subdwarf fitting (e.g., Reid 1997, 1998; COO). Note 



that KI03 derived cluster distances from eyeball fitting on 
(B - y, y), did not correct for a bias due to binaries (for both 
the cluster and subdwarf samples), and also did not apply 
the Lutz-Kelker corrections in their fitting procedure. We 
assumed lS.Ry = ±0.2 for the error in the reddening laws, 
and took a 20% error in E{B - V) and in the cluster age 
(Alogf (Gyr) « 0.1). We estimated the error from [a/Fe] by 
examining the case of [Q:/Fe]= 0.4 at [Fe/H] < 1.0 (§|3j. 

Systematic errors in the comparison of model colors are 
listed in Table |6] They are the average values estimated 
from all of our sample clusters. To avoid a large discrep- 
ancy in the models at the lower MS for M67 and NGC 6791 
(see below), we took the average in the color difference over 
4500 < reff(K) < 6500 for each cluster In addition to the in- 
dividual error sources described above, we also included an 
error from the fitting residual. Even after shifting isochrones 
in colors, we often found that they do not perfectly match the 
shape of the fiducial sequence. We computed this as an er- 
ror in the average color difference between the model and the 
cluster fiducial sequence. We discuss the effect of the globular 
cluster [Fe/H] scale in more detail in the next section. The to- 
tal error is the quadratic sum of all of the error contributions, 
which is dominated by the error in the cluster distance. 

As shown in Figures IgfflTl there is no compelling evidence 
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[Fe/H] = -2.4 [Fe/H] = -1.6 [Fe/H] = -1.2,-0.8 [Fe/H] = +0.0, + 0.4 




0.3 0.6 0.3 0.6 0.3 0.6 0.9 0.4 0.8 1.2 



(g - r)o (g - r)o (g - r)o (g - 0^ 

Fig. 5. — Reddening corrected (g — r,M,-) CMDs for sample clusters. Open and closed circles are fiducial points for each cluster. Left three panels: Solid lines 
are YREC+MARCS models at the age of 12.6 Gyr with [Fe/H] = -2.4,-1.6,-1.2,-0.8, as indicated on top of each panel. Distances from the Hipparcos-haseA 
subdwarf fitting are assumed, with canonical cluster reddening values (see text). The fiducial points for M71 are from Clem et al., after transforming u' g'r'i'z' 
into ugriz using transformation equations in Tucker et al. Right panel: Solid lines are models with [Fe/H] = +0.0 (4.0 Gyr), +0.3 (10.0 Gyr). Bull's-eyes mark 
the approximate position of the Sun, to guide eyes on the metallicity sensitivity of colors and magnitudes. 
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[Fe/H] = -2.4 [Fe/H] = -1.6 [Fe/H] = -1.2,-0.8 [Fe/H] = +0.0, + 0.4 




0.0 0.4 0.8 1.2 0.4 0.8 1.2 0.4 0.8 1.2 1.6 0.5 1.0 1.5 2.0 

(g - z)o (g - z)o (g - z)o (g - z)o 

Fig. 7.— Same as in Fig.|5] but (g-z,Mr). 
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TABLE 6 

Systematic Errors in the Comparison of Model Colors 



Source of EiTor 


A Quantity" 


A(g-r) 




Mg-z) 


A(u-g) 


Fitting Residual 




±0.003 


±0.004 


±0.006 


±0.008 


Phot. Calibration 


l%-2% 


±0.007 


±0.007 


±0.008 


±0.011 


(m-M)o 


±0.15 


±0.020 


±0.029 


±0.031 


±0.017 


Rv 


±0.2 


±0.002 


±0.002 


±0.001 


±0.006 


E(B-V) 


20% 


±0.006 


±0.012 


±0.019 


±0.008 


logt 


20% 


±0.003 


±0.005 


±0.009 


±0.008 


[a/¥ef 


±0.1 


±0.002 


±0.003 


±0.003 


±0.002 


[Ve/Hf 


±0.1 


±0.007 


±0.009 


±0.009 


±0.014 


Total 




±0.024 


±0.034 


±0.040 


±0.029 



Note. — Values are shown for the average color difference over 4500 < 7Ifr(K) < 6500. 
^ Approximate size of the enors. See §[2]for details.*^ Open cluster data aie not included in 
the eiTor estimation. 
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Fig. 9. — Color differences m g-r between YREC+MARCS models and 
fiducial sequences. The KI03 [Fe/H] values are assumed with the Hipparcos- 
based subdwarf fitting distances for globular clusters. Individual cluster cases 
are shown in such a way that the most metal-poor cluster is shown at the top, 
while the most metal-rich open cluster is shown at the bottom, with an offset 
A(g-r) = 0.10 between them. Dotted lines indicate zero differences between 
model and the data. For each cluster, filled and open circles exhibit the size 
of a systematic error from the photometric calibration (see text). A grey strip 
represents a ±1(t range of a total systematic eiTor in the compaiison. The 
result for M7 1 is based on Clem et al. fiducial sequence (after a photometric 
transformation to ugriz). 



that model colors differ from the observed MS for globular 
clusters in g-r, g-i, and g-z- On the other hand, there is a 
strong departure of the models at the lower MS for two open 
clusters {T^ff < 4500 K). The models predict bluer colors than 
the observed ones in all color indices. In fact, the problem 
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Fig. 10. — Same as in Fig.|9] but model comparisons in g — i. 



has long been noted in other photometric systems, such as 
the Johnson-Cousins system (e.g., S andquist.2004.) . for many 
contemporary low-mass stellar models. We present empirical 
corrections for the models in § 14.41 

There is an 0.1 mag offset in u-g for all clusters as 
shown in Figure [12] with no apparent trend in metallicity. 
Since the other color indices based on the g bandpass exhibit 
good agreement with the data, the u-g color offset suggests 
a problem in the u bandpass. The ^ 5% color difference is 
larger than expected from the photometric calibration error 
in SDSS. Instead, the difference is most likely due to miss- 
ing opacities in model atmospheres at short wavelengths. The 
difference is in the sense that models overestimate the flux in 
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Fig. 11. — Same as in Fig.|9] but model comparisons in g—z- 
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Fig. 12. — Same as in Fig.|9] but model comparisons in u—g. 



the u bandpass. The comparison in u-g for M71 is worse 
than those for the other clusters in our sample, suggesting a 
problem in the color transformation for this highly reddened 
cluster. 

4.2. Metallicity Scale for Globular Clusters 

Figure [13] shows the color differences between 
YREC+MARCS models and fiducial sequences, as a 
function of cluster metallicity, on the KI03 scale. Each 
open circle represents an average color difference for indi- 
vidual clusters over 4500 < 7'eff(K) < 6500. The grey strip 
displays a ±1(T error in Table |6l excluding the error from 
[Fe/H]. As seen in the figure, model colors are generally 
bluer than fiducial sequences for globular clusters. The 
average color differences are A(g-r) = -1-0.008 ± 0.01 1, 
A(^-/) =-1-0.022 ±0.017, A(g-z) = -1-0.021 ± 0.018, and 
A{u-g) = -1-0.059 ± 0.009 for the five globular clusters in our 
sample. 

Figure [14] shows a model comparison when the CG97 
[Fe/H] values are used. In comparison to the KI03 case 
(Fig. [T3T l. color differences between our models and fiducial 
sequences are reduced by ^ 0.02 mag: A{g- r) = -0.006 ± 
0.011, A(g-/)=-i-0.004±0.017, A(g-z)=-i-0.002±0.018, 
and A{u-g) = -1-0.031 ± 0.010. As discussed in KI03, their 
[Fe/H] values are about 0.2 dex lower than those from CG97, 
which make the isochrones redder. 

We took the difference from the two [Fe/H] scales as an 
effective 2a error, which is listed in Table [6] However, the 
metallicity scale has a relatively small impact on the model 



comparison, since the total error budget is still dominated by 
the error in the subdwarf-fitting distances. In the above anal- 
yses, we did not take a random error in [Fe/H] (cr^pe/Hj) ~ 
±0.01-0.02) into account, as it is much smaller than the sys- 
tematic error (cr^jFe/Hj) ~ ±0.1). We note that both CG97 and 
KI03 were based on spectroscopic measurements of red gi- 
ants, which could be different from those for dwarfs. How- 
ever, the difference m ay not be larger than ~ 0.1 dex (e.g., 
lKornetalJl2006ll2007h . 

4.3. Tests of Other Models 

Figure [15] shows a color comparison between cluster fidu- 
cial sequenc e s and Padova+ATLAS9 isochrones. Because 
iGirardi et alj (|2004|) did not include a-element enhanced 
models, we used their models in the comparison at the fol- 
lowing [Fe/H] values, in order to bracket the observed [Fe/H] 
abundance of a cluster (e.g., ISalaris et al.l 119931: IKim et aP 
I2OO2I) : [Fe/H] = -2.3 and -1.7 for M15 and M92, [Fe/H] = 
-1.7 and -1.3 for M3 and M13, [Fe/H] = -1.3 and -0.7 for 
M5, and [Fe/H] = -0.7 and -0.4 for M71. For the open 
clusters, we used their models at [Fe/H] = -l-O.O for M67, 
and [Fe/H] = -1-0.37 for NGC 6791. We appHed AB correc- 
tions for SDSS magnitudes (§ [3]). As shown in Figure [15] 
Padova+ATLAS9 models describe the observed cluster data 
reasonably well. However, the error bars (grey strip) are larger 
than those for the YREC+MARCS cases (Figs. [H] and [Ttl). 
because of the wider range of [Fe/H] values adopted in the 
model comparison. 

Similarly, Figure [16] shows a comparison between 
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Fig. 13. — Average color difference between the models and fiducial se- 
quences in ugriz. Open circles represent individual clusters, showing an av- 
erage color difference computed at 4500 < T^(((K) < 6500. For globular 
clusters the KI03 [Fe/H] values are used. Error bars and grey strips represent 
±1(T range of a total systematic error. 

DSEPh-PHOENIX models and fiducial sequences. Since their 
models are available over a wide range of [Fe/H] and [a/Fe], 
isochrones were interpolated at the cluster's [Fe/H] on the 
K103 scale. Although model colors and magnitudes from 
DSEPh-PHOENIX already included AB corrections, 0.01 mag 
was subtracted from their g-band magnitudes, to be consis- 
tent with our adopted AB corrections for SDSS. Although 
DSEPh-PHOENIX models are in reasonably good agreement 
with the cluster data, they show a large difference in g- 
z and u- g from the data. The above tests suggest that 
YREC+MARCS models are generally in better agreement 
with observed cluster data than these models. 

4.4. Empirical Corrections on Color-T^tf Relations at 
[Fe/H]> 

In the above two sections, we restricted our model compari- 
son to Teff > 4500 K because the models are in poor agreement 
with the lower MS at solar and super-solar metallicity (Fig.|9]- 
[T2li. Both the Padova+ATLAS9 and DSEP+PHOENIX 
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Fig. 14. — Same as in Fig. [13] but when the CG97 [Fe/H] values are used 
in the model comparison for globular clusters. 

isochrones show similar discrepancies in colors for open clus- 
ters. Since CMDs for globular clusters do not extend far be- 
low Teff ~ 5000 K, similar problems could not be found at 
the lower MS. However, the comparison with M71 suggests 
that model colors and magnitudes are probably accurate at 
[Fe/H] < -0.8. 

Figure [T7] shows color differences for individual cluster 
members in M67. These stars were selected using photo- 
metric filtering of background stars and cluster binaries on 
(g-r,r), {g-i,r), and (g-z,r), independent from theoretical 
models (see Paper 1). We used a solar-metallicity isochrone 
at an age of 3.5 Gyr For our calibrati on exercise in the 
Johnson-Cousins-2MA SS filter system (iPinsonneault et al] 
120041: lAn etal.ll2"007al) . we used photometry of the Hyades 
and Pleiades open clusters to define empirical corrections on 
the color- Teff relations. Similarly, we could define empirical 
corrections on the model colors in ugriz based on M67. 

The sohd fines in Figure [T7] are polynomial fits to the 
points: 



A(^-r) = 7.613- 



■ 19.91061+17.27961^-4.983611 



(3) 
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Fig. 15.— Same as in Fig. [HI but for Padova+ATLAS9 models. For 
each globular cluster, isochrones with two different [Fe/H] values are used in 
the model compaiison to bracket the observed cluster abundance (see text), 
because Padova+ATLAS9 models do not include a-element enhancement. 

A(g-0 = 6.541- 16.31361-1- 13.4616'2-3.6756'^ (4) 

A(g - z) = 3 .755 - 7.59561 -i- 4.5356(2- 0.6726l^ (5) 

A(m - = -7 .923 -I- 24. 1 3 1 61 - 23 .7279^ +7.6279^, (6) 

where 9 = T^n/ 5040 K. These corrections are in the sense that 
the above values should be added to the model colors, and 
they are strictly valid at 4000 < 7;ff(K) < 6000. 

It should be noted that the underlying assumptions used in 
these cor rections are somewhat different from our previous 
work. In lPinsonneault et all (I2003L |2004|) . the main cause of 
the residual difference between the model and the observed 
MS was attributed to the model color- Teff relations, since the 
mass-luminosity-Teff relations were found valid over 4500 < 
7eff(K) < 6000, where the color corrections were defined. On 
the other hand, the discrepancy between the lower MS and 
the models at T^a < 4500 K in Figure [17] could also be due 
to our limited knowledge of the opacities and/or the equations 
of state in the models at low temperatures, in addition to a 
problem in the model atmospheres. 

We tested the above corrections using CMDs for 
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Fig. 1 6 .— Same as in Fig.[73] but for DSEP-hPHOENIX models. TheKI03 
[Fe/H] values are used for globular clusters. 

NGC 6791, as shown in Figure [18] Dotted lines are 
YREC+MARCS isochrones without corrections, and solid 
lines are those with M67-based color-Teff corrections. Al- 
though [Fe/H] -dependent corrections can be derived using 
NGC 6791, empirical corrections based on the solar metal- 
licity cluster result in a reasonably good agreement with the 
observed MS for the super-solar metallicity cluster The cali- 
brated isochrones appear to deviate from the observed MS at 
the faint end, but the location of the observed MS becomes 
uncertain as well, due to the increased scatter Because of a 
paucity of clusters at -1 < [Fe /H] < in the sample, the lower 
metallicity limit of the color corrections is not well defined. 

5. DISTANCES AND AGES OF GLOBULAR CLUSTERS 

The ages of Galactic globular clusters have been ex- 
tensively discussed in the literature, because of their im- 
portance in setting a lower limit on the age of the uni- 
verse, and their im plications for t he formation of the 
Galact ic halo (e.g., Sandage 1970; Searle & Zinn 19781 
Chaboveret al.1 11998; Rosenberg et al. 1999; YandenBer'3 
20001; ISalaris & Weissll2002l; iDe Angeli et al.ll2005h . Recent: 
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Fig. 17. — Empirical corrections on theoretical color- Teff relations. Points 
represent color differences of individual stars in M67 with respect to a solar- 
metallicity, 3.5 Gyr old YREC+MARCS isochrone. These stars are those 
selected from a photometric filtering as single MS members based on (g — 
r,r), (g — i,r), and {g-z,r). Solid lines are empirical color corrections to the 
models, which are the polynomial fits to the points in each color index. 

age estimates for the Galactic globular clusters are typi- 
cally found between ~ 10 Gyr and ~ 15 Gyr based on var- 
ious techniques, including MSTO fitting (e.g., iGratton et al.l 
[19971 [Chaboyer et al."1998'; 'VandenBerg'2000'), termination 
of the w hite dwarf cooling sequence (e.g. Renzini et al. 1996; 
IZoccali et a l. 2001; Hansen et al.ll2002l) . luminosity function 
(e.g.l Jimenez & Padoan 1998|), and the mass-age relati on for 
an eclipsing binary system IChabover & Kraussll2002l) . Al- 
though relative age estimates for globular clusters have been 
well studied (e.g.. Stetson et al. 1996; Rosenberg et al. 199^, 
estimating absolute ages is a more difficult problem, primar- 
ily b ecause of the uncertainty in distances to globular clusters 
(e.g.. lVandenberg et al.lll996l) . 

In this section, we estimate distances to our sample globu- 
lar clusters using the MS -fitting technique without employing 
empirical color corrections. We then use the derived distance 
to estimate cluster ages from the MSTO. We test the valid- 
ity of deriving these quantities from the SDSS photometric 



database by comparing our values with those in the literature 
using other methods. In particular, we show that our distances 
using YRECh-MARCS models are in good agreement with the 
Hipparcos-hsised subdwarf-fitting distance estimates. 

5.1. Distances to Globular Clusters 

We obtained distances to globular clusters from CMDs with 
three color indices (g-r, g-i, and g-z). We did not use 
u-g because of the steep MS slope in the {u-g,r) plane 
(Fig. [8]l. The utility of the mu lti-color fitting technique has 
been demonstrated in I An et aTl (2007ajb). With a wide wave- 
length baseline, stellar properties are not only accurately con- 
strained, but systematic errors can be estimated. 

Table |7] lists the derived distance moduli for our sample 
clusters, using fiducial sequences on (g - r, r), {g - i, r), and 
(g-z,r). We restricted our fits to AM^ k, 1 mag below the 
MSTO, to minimize any age dependency. For the globular 
clusters, we list four different sets of distance estimates, on 
two different metallicity scales, and with two different sets 
of photometric zero points (§ |4]l. We computed the size of 
the error in the average distance modulus from the propaga- 
tion of a fitting error (shape mismatch) or from a dispersion 
of distance moduli from three CMDs, whichever is larger (see 
below for other systematic errors). We found, however, that 
the latter is always larger than the former by approximately 
a factor of four The last column lists a standard deviation 
of the distance moduli from three CMDs, indicating that dis- 
tances from individual CMDs are consistent with each other 
to better than ^ 5%. In Table |7] we also included results for 
open clusters (M67 and NGC6791) usi ng p urely theoretical 
models without empirical corrections (§ I4.4| | at r < 16.4 mag 
and r < 21 .0 mag for M67 and NGC 6791, respectively. 

Table |8] lists systematic errors in distance estimation. We 
used the same list and size of systematic error sources as in 
Table |6] An "internal" error was taken from the error in the 
average distance modulus in Table |7] Although the adopted 
size of the photometric calibration errors varies from one clus- 
ter to the other, the small size of the errors (e.g., M13) may 
reflect an accidental agreement of photometry from two SDSS 
runs. The error from [Fe/H] increases at higher metallicities, 
because the metallicity sensitivity of colors and magnitudes 
increases steeply with [Fe/H] . A quadrature sum of all of these 
systematic errors is generally larger than the internal system- 
atic error. We note that no particular source of systematic 
errors dominate the total error budget in our distance determi- 
nation. Average distance moduli for individual clusters and 
their total errors are listed in the last two columns of Table |2] 
on the CG97 and K103 metallicity scales, respectively. 

Figure [19] compares our distance moduli for globular clus- 
ters with those based on Hipparcos subdwarf fitting. Open 
circles are those estimated using the CG97 [Fe/H] values, 
while closed circles are those estimated using the K103 [Fe/H] 
values. Figure l20l shows these differences, as a function of 
metallicity, on the two different [Fe/H] scales. Error bars are 
the quadrature sum of the errors from both our values and 
those from the subdwarf-fitting studies. The average differ- 
ence between these two studies is A(m-M)o = -0.02 ±0.07 
on the CG97 [Fe/H] scale, and A(m-M)o = -0.08 ±0.07 on 
the K103 scale. In both cases, our distances are consistent 
with subdwarf-fitting distances. Furthermore, the agreement 
between subdwarf-fitting and this study indicates that our er- 
ror estimate ((T(,„_m) ^ 0.03-0. 1 1 mag) for individual clusters 
is reasonable. In fact, the agreement between the two dis- 
tance determinations restates our earlier conclusion that the 
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Fig. 18. — CMDs for NGC 6791 . D otted lines are theoretical models with [Fe/H]= +0.4 at the age of 10 Gyr. SoHd Hnes are models with empirical color- Teff 
corrections based on M67 (see Fig. 117) . 



TABLE 7 
MS-FiTTiNG Distances 





GC 


Reference 






(m-M)o 






Cluster 


[Fe/H] 


SDSS Run" 


(8-r,r) 


(S->,r) 


{g-z,r) 


Average 


O" (m-M) 


M15 


K103 


2566 


15.061 ±0.004 


15.061 ±0.007 


14.930 ±0.039 


15.060 ±0.053 


0.092 




K103 


1739 


15.150±0.016 


15.142±0.011 


15.011 ±0.048 


15.140±0.053 


0.091 




CG97 


2566 


15.101 ±0.007 


15.094±0.008 


14.958 ±0.043 


15.096 ±0.056 


0.098 




CG97 


1739 


15.190±0.019 


15.176±0.014 


15.038 ±0.052 


15.175±0.056 


0.097 


M92 


K103 


4682 


14.619±0.019 


14.577 ±0.013 


14.605 ±0.016 


14.599 ±0.012 


0.021 




K103 


5327 


14.745 ±0.013 


14.714 ±0.006 


14.549 ±0.007 


14.655 ±0.062 


0.107 




CG97 


4682 


14.654 ±0.017 


14.604 ±0.012 


14.628 ±0.011 


14.624 ±0.015 


0.026 




CG97 


5327 


14.780 ±0.014 


14.741 ±0.007 


14.571 ±0.006 


14.656 ±0.070 


0.122 


M13 


K103 


3225 


14.357 ±0.007 


14.346 ±0.003 


14.326 ±0.004 


14.341 ±0.009 


0.016 




K103 


3226 


14.393 ± 0.006 


14.326 ±0.005 


14.3 16 ±0.005 


14.340 ± 0.024 


0.042 




CG97 


3225 


14.441 ±0.010 


14.417 ±0.006 


14.392 ±0.008 


14.414 ±0.014 


0.025 




CG97 


3226 


14.477 ±0.010 


14.398 ±0.006 


14.382 ±0.005 


14.400 ±0.032 


0.056 


M3 


K103 


4646 


14.822 ±0.011 


14.884 ±0.003 


14.851 ±0.028 


14.879 ±0.026 


0.045 




K103 


4649 


14.995 ±0.006 


14.981 ±0.004 


14.938 ±0.032 


14.985 ±0.020 


0.034 




CG97 


4646 


14.921 ±0.009 


14.971 ±0.005 


14.939 ±0.034 


14.959 ±0.018 


0.031 




CG97 


4649 


15.093 ±0.010 


15.068 ±0.009 


15.025 ±0.039 


15.078 ±0.023 


0.039 


M5 


K103 


1458 


14.329 ±0.006 


14.239 ±0.002 


14.294 ±0.007 


14.251 ±0.037 


0.063 




K103 


2327 


14.353 ±0.009 


14.287 ±0.005 


14.247 ± 0.003 


14.265 ±0.038 


0.066 




CG97 


1458 


14.426 ±0.010 


14.326 ±0.002 


14.376±0.012 


14.331 ±0.043 


0.074 




CG97 


2327 


14.450 ±0.012 


14.373 ±0.007 


14.328 ±0.006 


14.360 ± 0.039 


0.068 


M71 


KI03 




12.846 ±0.022 


12.786 ±0.031 


12.898 ±0.045 


12.836 ±0.033 


0.057 




CG97 




12.946 ±0.028 


12.875 ±0.034 


12.980 ±0.047 


12.928±0.031 


0.054 


M67 




5935 


9.633 ±0.014 


9.584 ±0.014 


9.572 ±0.011 


9.592 ±0.019 


0.033 






5972 


9.680 ±0.021 


9.578 ±0.013 


9.601 ±0.013 


9.604 ±0.033 


0.057 






6004 


9.637±0.114 


9.588 ±0.014 


9.620 ±0.014 


9.615±0.014 


0.025 


NGC 6791 




5416 


13.005 ±0.021 


12.944±0.015 


12.932 ±0.014 


12.951 ±0.024 


0.041 






5403 


12.930 ±0.016 


12.954±0.014 


13.013±0.019 


12.960 ±0.025 


0.043 






6177 


13.000 ±0.019 


12.962 ±0.014 


12.971 ±0.015 


12.974 ±0.012 


0.020 



SDSS runs selected for the local photometric zero points. 
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TABLE 8 

Systematic Errors in Distance Moduli for Individual Clusters 



Source of 
Error 




















A Quantity 


M15 


My2 


M13 


M3 


M5 


M71 


M67 


NGC 6791 


Internal 




±0.053 


±0.012 


±0.009 


±0.026 


±0.037 


±0.033 


±0.019 


±0.024 


Phot. Calibration 


l%-2% 


±0.040 


±0.028 


±0.000 


±0.053 


±0.007 


±0.026 


±0.012 


±0.012 


Rv 


±0.2 


±0.009 


±0.002 


±0.001 


±0.002 


±0.008 


±0.040 


±0.003 


±0.006 


E{B-V) 


20% 


±0.089 


±0.013 


±0.043 


±0.008 


±0.019 


±0.055 


±0.010 


±0.032 


\ogt 


15% 


±0.014 


±0.009 


±0.007 


±0.014 


±0.039 


±0.003 


±0.007 


±0.008 


[a/Fe] 


±0.1 


±0.008 


±0.005 


±0.017 


±0.011 


±0.020 


±0.015 


±0.000 


±0.000 


[Fe/H] 


±0.1" 


±0.018 


±0.013 


±0.037 


±0.040 


±0.040 


±0.046 


±0.012 


±0.027 


Total 




±0.114 


±0.037 


±0.060 


±0.074 


±0.073 


±0.093 


±0.028 


±0.051 



^ Error in the metallicity scale between CG97 and K103 for globular clusters. The (7[pe^H] = iO.Ol and cr[Fe/H] = ±0.03 were assumed for 
M67 and NGC 6791, respectively. 



Galactic Star Clusters in SDSS II. 



19 



15.0 



m 

m 
•I— I 



14.0 - 



S 

I 



13.0 

0.4 
0.0 
-0.4 



"1 ' ' ' ' r 

O CG97 [Fe/H] 
• KI03 [Fe/H] 



H ^ — \ — h 



H — I — \ — h 



H h 



_L 



_L 



_L 



13.0 14.0 15.0 

(m - M)o (Subdwarf Fitting) 

Fig. 19. — Comparison of distance moduli for globular clusters from the 
Hipparcos-hased subdwarf fitting (see § |2) and this study. Open circles are 
those estimated from the K103 [Fe/H] values, and the closed circles are those 
from the CG97 [Fe/H] values. 
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Fig. 20. — Same as in Fig. [19] but as a function of metallicity. Top: 
Comparison of distance moduH when the CG97 [Fe/H] values are adopted 
in the MS fitting. Bottom: Comparison when the KI03 [Fe/H] values are 
adopted. 



YREC+MARCS models are in good agreement with cluster 
fiducial sequences. 

The distance modulus of M71 from lGrundahl et all (l2002h 
is 0.33 ±0.17 mag smaller than the value in iReidld 19981) 
although these studies assumed the same E(B-V) « 0.28. 
Our best-fitting distance modulus of this cluster, (m-M)o = 
12.86 ±0.08 (KI0 3) or 12.96 ±0 .08 (CG97), are in better 
agreement with the lGrundahl et alj value. Our cluster distance 
modu lus is also in agreemen t with (m-M)o = 12.76±0.18 
from lGeffer t & Maintz (2000). In addition, the average value 
from the Harris compilation becomes (m-M)o = 12.92 when 
we use E(B-V) = 0.28, rather than his adopted E(B-V) w 
0.32. 

COO paid particular attention to M92, not only because of 
its importance as one of the most metal-poor Galactic clus- 
ters, but also because of the large discrepancies in its re- 
ported distance in the literature. Two factors make measur- 
ing the cluster's distance difficult. First, there are only a few 
metal -poor subdwarfs with [Fe/H] < -2 in the Hipparcos 
catalog. Secondly, the two independent photometric studies 
(Heasley & Christian 1986; Stetson & Harris 1988) that pro- 
vided BV photometry have 0.03 mag differences in colors. As 
shown in Table |2] with independent photometric data and a 
wide wavelength coverage in griz, we found that our distance 
to M92 is in good agreement with the COO distance, indepen- 
dent of the metallicity scale used. 

Here w e neglected the effects of unresolved binaries. In 
lAn et al.l (l2007h) we have performed extensive simulations of 
unresolved binaries in clusters and inspected their influence 
on the estimation of a MS-fitting distance. After applying 
photometric filtering, as was done in the present paper for our 
cluster sample, we found that unresolved binaries can make a 
MS look brighter by ~ . 007 mag for a 40% binary fraction ' . 
The observed bina ry fraction in globular clusters is typi cally 
less than 20% (e.g.. lSomma et al.l2007l:rDavis et alj20i08 , and 
references therein), which makes the influence of unresolved 
binaries even smaller. 

5.2. Cluster Ages 

We have derived cluster ages for our sample using MSTO 
stars, based on the MS-fitting distances to clusters described 
above. Figures I2TH251 displav {g-r,r), (g-i,r), and {g-z,r) 
CMDs, respectively, for five globular clusters in our sample, 
with best-fitting stellar isochrones (solid lines, see below). We 
fit directly to the CMDs, instead of using the cluster fiducial 
sequences, because fiducial sequences can be inaccurate near 
the MSTO (where curvatures are largest on CMDs). 

The first panel in Figures |2TI - |25] exhibits residuals on (g- 
r, r) with respect to a best-fitting model. The points in the re- 
maining three panels are those after applying cuts based on 
the X, sharp, and separation indices, as well as the photomet- 
ric filtering (see Paper I for more details). The photometric 
filtering routine statistically eliminates cluster non-members 
and unresolved binaries, independent of the stellar isochrones. 
However, we note that cuts based on the separation index al- 
ready make a cluster fiducial sequence sufficiently naiTow for 
well-populated globular clusters such as those in our sample. 

To find a best-fitting age on each CMD, we minimized 



Xtot ■ 



(Xi X,„) 



(7) 



" Binary fraction is defined as the number of binaries divided by the total 
number of systems. 
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Fig. 2 1 . — CMDs for M15 with YREC+MARCS models at the best-fitting distance and age (solid line). Models without empirical color corrections are used. 
Dotted, horizontal lines represent a magnitude range where the isochrone fit was performed to estimate the cluster age. Points are stars that were remained after 
the photometric filtering (see text). The leftmost panel shows the fitting residuals in (g — r,r). Histograms are the number distribution of the residuals, with an 
arbitrary scale on the vertical axis. 
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Fig. 22.— Same as in Fig.ET] but for M92. 
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Fig. 23.— Same as in Fig.[2T] but for M13. 
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TABLE 9 
MSTO Ages of Globular 
Clusters 



Cluster 


Age 


(Gyr) 




Lhe/HJcG97 


M15 


13.9±2.5 


12.8±2.3 


M92 


14.4±0.9 


13.5 ±0.9 


M13 


14.3±1.I 


12.3 ±1.0 


M3 


13.3±1.4 


11.4±1.1 


M5 


12.2±1.3 


10.5 ±1.2 



at the best-fitting distance derived simultaneously using MS 
stars. Here, X, and ax,i are a color and its error for the 
point, respectively. The CTr,,- is an error in r-band magnitude, 
weighted by the inverse slope of the isochrone 7. The X,„ is 
a model color at the star's r-band magnitude. We used stars 
that are within ^ 0.5 mag in r on the bright side of the bluest 
point of the MSTO and 1.0 mag on the faint side. The 
bright magnitude limit was set by the fact that our evolution- 
ary model computations do not reach as far beyond the hydro- 
gen exhaustion in the core. On the red side of subgiant branch 
(SGB), assumptions on the mixing length and/or the model 
color- Teff relations become uncertain as well. The faint end 
was set at this level, since the MS is insensitive to the age be- 
low this limit. Dotted horizontal lines in Figures |2T]-|25] show 
these fitting ranges. 

Solid lines in Figures |2T]-|25]represent our best-fitting mod- 
els for individual CMDs, using KI03 [Fe/H] values. The g-r 
residuals in the first panels show that fits are excellent from 
the lower MS to the MSTO. Histograms are also shown to 
visualize a peaked distribution of color residuals from the 
model. Reduced values were between ~ 0.9 and ~ 1.2 
for approximately 2000 stars on each CMD. 



Our best-fitting ages are tabulated in Table |9] The errors 
are the quadrature sums of various systematic errors in the es- 
timation, as listed in Table [TOl The same sizes of systematic 
errors were adopted as in Tables |6] and |8] The first row in 
Table [TO] shows the scatter in the age estimation from three 
CMDs (g-r, g-i, and g-z), which describes an internal sys- 
tematic error of ~ 0.5 Gyr. Since the distance and age esti- 
mates are correlated with each other, we evaluated the effects 
of various systematic errors, while simultaneously solving for 
both quantities. In Table [TOl we separately tabulated the con- 
tribution from an internal error component in distance (see 
Table [U. We took the difference in age estimates from the 
CG97 and KI03 [Fe/H] scales as a 2a error Our error esti- 
mates indicate that we can obtain cluster ages with an accu- 
racy of ~ 10%-15% using multiple color indices in SDSS. 

The accuracy of our age estimates can be further checked 
through a comparison with previous studies. Figurel26l shows 
compari sons of our a ge estimates wit h those in two recent 
studies: 'VandenBere (2000, top) and ISalaris & WeissI (I2002L 
bottom). VandenBerg (2000) estimated cluster ages from the 
luminosity difference between his theoretical HB models and 
two points on the cluster's MS and SGB, which are 0.05 mag 
redder in B-V than the bluest point of the MSTO. Be- 
cause the cluster [Fe/H] values adopted in l VandenBerd (l2000h 
were generally found to lie between the CG97 or KI03 val- 
ues,^" we used the average of our age estimates at the CG97 
and KI03 [Fe/H] values in the comparison. On the other 
hand, ISalaris & \\feissl (l2002h derived ages for M3 and Ml 5 
from the V magnitude difference between Zero- Age HB and 
MSTO, but employed relative age estimates for other clusters 
in our sample from the color difference in B - V or V - /c be- 

^° I VandenBerd tSOOOh adopted [Fe/H]= -2.30 for M15 and M92, [Fe/H]= 
-1.60 for M3 and M13, and [Fe/H]= -1.40 for M5. 
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TABLE 10 

Systematic Errors in the Ages of Globular Clusters 



Source of 
EiTor 








""log I 






A Quantity 


M15 


Mn 


M13 


M3 


M5 


Internal 




±0.024 


±0.014 


±0.009 


±0.012 


±0.027 


Phot. Calibration 


l%-2% 


±0.022 


±0.015 


±0.006 


±0.024 


±0.004 


(m-M)o 


Internal 


±0.003 


±0.003 


±0.001 


±0.003 


±0.002 


Rv 


±0.2 


±0.009 


±0.003 


±0.002 


±0.000 


±0.002 


E(B-V) 


20% 


±0.061 


±0.012 


±0.019 


±0.007 


±0.019 


[Q/Fe] 


±0.1 


±0.010 


±0.011 


±0.023 


±0.029 


±0.028 


[Fe/H] 


±0.1 


±0.018 


±0.015 


±0.032 


±0.033 


±0.034 


Total 




±0.073 


±0.030 


±0.045 


±0.052 


±0.055 






(±16.8%) 


(±7.0%) 


(±10.4%) 


(±12.0%) 


(±12.7%) 
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Fig. 26. — Top: Comparison of globular cluster age estimates in 
l^ndenBera (2000) and those in this study. The latter are the average val- 
ues computed from the CG97 and K103 [Fe/H] values, while the former age 
estimates are based on those [Fe/H] values adopted in that study (see text). 
Bottom: Comparison between Salaris & Weiss. ( 2002) and this study. Both of 
these age estimates are based on the CG97 [Fe/H] values. 



tween the MS TO and the base of R GB. They presented ages 
from both the lZinn&\\^ digs! and CG97 [Fe/H] scales, 
but only those from the latter [Fe/H] scale are shown in Fig- 
ure |26] 

There is excellent agreement between our age estimates and 



those in I VandenBerd ('2000') (Fig. |26] |. even though they em- 
ployed independent photometry, photometric filters, stellar 
isochrones, and different methods of estimating cluster dis- 
tances and ages. In addition to the above two studies, COO de- 
rived ages for 9 globular cluster s by comparing MSTO lumi- 
nosities with model isochrones (IStraniero et al.lfT997l) . They 
obtained 14.8, 12.6, and 11.2 Gyrs for M92, M13, and M5, 
respectively, which are in agreement with our values. The 
outcome of the comparisons between ours and the preceding 
studies demonstrates that cluster ages derived with the SDSS 
ugriz system are reliable, and that our error estimates on age 
are reasonable. 

However we also noted above that theoretical er- 
rors (e.g. [V andenberg eLalJ 119961: IChabover et al] 119981: 
iKrauss & C habover 2003) were not included in our total error 
budg et. For example, the inclusion of the microscopic diffu- 
sion dKorn et al. 2006, 2007) can systematically reduce stellar 
ages by 10% (e.g..lChabover et alj200ll;IVandenBerg et alj 
l2002t iMichaud et al Il2004 . Nevertheless, our fitting results 
(Figs. |2T1 - |25] | are encouraging, showing excellent fits in all 
three of the color indices. This suggests that our mixing length 
calibration based on the Sun is valid for stars in globular clus- 
ters, and that theoretical color- Teff relations are also good for 
MSTO stars, as we found for MS stars in §|4l 

Cluster ages have traditionally been estimated from the 
absolute V magnitude of the bluest point on MSTO (e.g., 
Ay,^ method). This was because the luminosity of the 
MSTO was considered to be the most robustly predicted quan- 
tity in stellar evolution ary models (e.g., B olte & Hoganlll995l : 
IChaboyer et alJ[T996bh . Our approach of fitting isochrones to 
the whole MSTO has been cautioned against in the literature 
(e.g., Renzini 1991), primarily because of limits on knowl- 
edge of the mixing length parameter and the color- Teff rela- 
tions used in the models. However, the method has 
its own limitations, in particular the observational difficulty 
of pinning down the bluest MSTO point. Because the MS 
is almost vertical in that region, the magnitude of the MSTO 
cannot be better determined than AV ^0.1 mag, which then 
produces an error of ^ 1.5 Gyr in age. Subsequent stud- 
ies modified the AV™ technique in such a way that slightly 
redde r points than the MSTO are used as an age indicator 
(e.g.. IChabover et all [19961 IVande nBerg 2000). However, 
this method implicitly assumes that the shape of MSTO in 
the models is correct, as we have in this paper 

Figure l27l shows MSTO ages for our sample globular clus- 
ters as a function of metallicity. The open circles are those 
estimated from the CG97 [Fe/H] values, while the filled trian- 
gles are those from the KI03 values. Error bars represent the 
total error in age, as listed in Table |9l With constraints on the 
distance, we can inspect how stellar ages can be affected by 
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Fig. 27. — Turn-off ages of globular clusters in this study as a function of 
metallicity. Filled triangles and open circles represent the cases of adopting 
the KI03 or CG97 [Fe/H] values, respectively. 

secondary factors, such as the cluster [Fe/H] scale. When the 
KI03 [Fe/H] values are adopted, the average age of our sam- 
ple clusters is 13.8 ±0.5 Gyr Some of the clusters are even 
older than the inferred age of the universe based on the Wilkin- 
son Microwave Ani s otwpy Probe ( WMAP) data (« 13.7 Gyr, 
ISpergel et al.1 120031: iDunkley et al.l l2008), although they are 
still consistent with each other within the errors. However, the 
average MSTO age becomes 12.2 ±0.5 Gyr when the CG97 
values are used, reducing tension between stellar and WMAP- 
based ages. In addition, the age estimates based on the CG97 
scale exhibit a steeper age-metallicity relation. Although only 
five clusters are included in the analysis. Figure [27] demon- 
strates that a different cluster [Fe/H] scale could lead to a dif- 
ferent chemical enrichment scenario for the Galactic globular 
cluster system. 

Finally, theoretical colors of the MSTO are presented in 
Figure |28] for three color indices, to aid the interpretation of 
CMDs in SDSS, as well as for future imaging surveys such 
as the Panoramic Su rvey Telescope & Rapid Response Sys- 
tem (Pan-STARRS; Kaiser et a l. 2002) and the Large Synop- 
tic Survey Telescope (LSST; Jvezic et al. 2008b), which will 
use similar ugriz filter system. No empirical corrections were 
applied to the theoretical colors. 

6. SUMMARY AND DISCUSSION 

We generated a set of stellar isochrones in the SDSS ugriz 
system and performed an extensive test of these models for 
MS stars. We found that the models match observed cluster 
sequences within the errors of the adopted cluster parameters, 
for CMDs using g-r, g-i, and ^-z as color indices. These 
models correctly predict stellar colors and magnitudes over a 
wide range of metallicity, given the size of errors in distance, 
reddening, metallicity, and photometric calibration. However, 
we found ~ 0.05-0.10 mag offsets in u-g, indicating a prob- 
lem in the model atmospheres for the u bandpass. More ac- 
curate tests wil l be possible when fut ure astrometric missions 
such as Gaia (iPerryman et al.l |200T1) provide improved dis- 
tance measurements for these clusters. 

On the other hand, we found a strong departure of the 
models from the open cluster observations of the lower MS 
stars. The problem with cool stars has long been noted in 
other photometric systems, such as in Johnson-Cousins band- 
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Fig. 28. — Theoretical MSTO colors in three color indices. 

passes, and is presumably due to our limited knowledge on 
the opacities and/or the equations of state used in the models. 
From the observation of M67, we proposed empirical cor- 
rections on the model color- Teff relations to alleviate a large 
distance error to the cluster from MS fitting. Our tests in- 
dicate that isochrones with these corrections provide an im- 
proved fit to MS of NGC 6791. Since colors and magnitudes 
of stars change rapidly with metallicity, accurate photometry 
of open clusters in ugriz will be useful to explore metallicity- 
dependent corrections on color- Teff relations in the future. 

We also tested models from Pado va-i-ATLAS9 
(iGirardi et al.l |2004) and DSEPh-PHOENIX (iDotter et all 
120081) using cluster fiducial sequences as for comparison 
with the YRECh-MARCS models - the models provide 
reasonable fits to the observations generally, but exhibit 
strong departures for the lower MS at solar and super-solar 
metallicities. However, Padova-i-ATLAS9 models do not 
include a-element enhancement, prohibiting more accurate 
comparisons with cluster observation. In addition to u-g, 
we also found a large offset in g-z for DSEPh-PHOENIX 
models. In general, our YRECh-MARCS models better match 
observed MSs than these models. 

We derived distances and MSTO ages for our sample 
of Galactic globular clusters using our theoretical stellar 
isochrones. We performed isochrone fits for three color in- 
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dices (g—r, g — i, and g-z) over a wide wavelength range 
to constrain these quantities and to assess the size of inter- 
nal systematic errors. Our distance estimates are consistent 
with those from Hipparcos-hased subdwarf fitting, with an 
error of (T(m-M) ^ 0.03-0.11 mag for individual clusters. We 
also found good agreement of our age estimates with those 
in the literature, with Cage ~ 10%-15%. Our results demon- 
strate the ability to derive distances and ages for stars ob- 
served in SDSS, as well a s in the future fr om imaging sur- 
yeys suc h as Pan-STARRS dKaiser e t al."2002) and with LSST 
dlvezic e t al. 2008b) . which will employ similar ugriz filter 
systems. 



We thank the referee for careful and detailed comments. 
D. A. thanks Grant Newsham for kindly proving his isochrone 
interpolation code. We thank Aaron Dotter and Brian 
Chaboyer for helpful discussions on the model comparison. 
D.A. gratefully acknowledges the support from the Presiden- 
tial Fellowship in the Ohio State University. F.D. is sup- 
ported by ANR (ANR-06-CIS6-009-01). T.C.B. is grateful 
for support in part by from grant AST 07-07776, and grants 
PHY 02-16783 and PHY 08-22648: Physics Frontiers Center 
/ Joint Institute for Nuclear Astrophysics (JINA), awarded by 
the U.S. National Science Foundation. Any opinion, finding, 
and conclusions or recommendations expressed in this mate- 



Abazajian, K., et al. 2003, AJ, 126, 2081 

Abazajian, K., et al. 2004, AJ, 128, 502 

Abazajian, K., et al. 2005, AJ, 129, 1755 

Adelman-McCarthy, J. K., et al. 2006, ApJS, 162, 38 

Adelman-McCarthy, J. K., et al. 2007, ApJS, 172, 634 

Adelman-McCarthy, J. K., et al. 2008, ApJS, 175, 297 

Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437. 879 

Alonso, A., Arribas, S., & Martmez-Roger, C. 1995, A&A, 297, 197 

Alonso, A., Arribas, S., & Martmez-Roger, C. 1996, A&AS, 117, 227 

An, D., Terndrup, D. M., & Pinsonneault, M. H. 2007a, ApJ, 671, 1640 

An, D., Terndrup, D. M., Pinsonneault, M. H., Paulson, D. B., Hanson, R. B., 

& Stauffer, J. R. 2007b, ApJ, 655, 233 
An, D., et al. 2008, ApJS, 179, 326 

Badnell, N. R., Bautista, M. A., Butler, K., Delahaye, F., Mendoza, C, 

Palmeri, P, Zeippen, C. J., & Seaton, M. J. 2005, MNRAS, 360, 458 
Bahcall, J. N., Pinsonneault, M. H., & Basu, S. 2001, ApJ, 555, 990 
Bessell, M. S. 1990, PASP, 102, 1181 
Bessell, M. S., Ca.stelli, F, & Plez. B. 1998, A&A, 333, 231 
Bolte, M., & Hogan, C. J. 1995, Nature, 376, 399 

de Bruijne, J. H. J., Hoogerwerf, R., & de Zeeuw, R T. 2001, A&A, 367, 1 1 1 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Carraro, G., Villanova, S., Demarque, P., McSwain, M. V., Piotto, G., & 

Bedin, L. R. 2006, ApJ, 643, 1151 
Carretta, E., & Gratton, R. G. 1997, A&AS. 121, 95 (CG97) 
Carretta, E., Gratton, R. G., Clementini, G., & Fusi Pecci, F 2000, ApJ, 533, 

215 (COO) 

Castelli, F, Gratton, R. G., & Kurucz, R. L. 1997, A&A, 318, 841 

Castelli, F, & Kurucz, R. L. 2003, in lAU Symp. 210, Modeling of Stellar 

Atmospheres, ed. N. E. Piskunov, W. W. Weiss, & D. F. Gray (San 

Francisco: ASP), poster A20 on enclosed CD-ROM I astro-ph/0405087 I 
Chaboyer, B., Demarque, P., Keman, P. J., Krauss, L. M., & Sarajedim, A. 

1996, MNRAS, 283, 683 
Chaboyer, B., Demarque, R, & Sarajedini, A. 1996, ApJ, 459, 558 
Chaboyer, B., Demarque, R, Keman, R J., & Kiauss, L. M. 1998, ApJ, 494, 

96 

Chaboyer, B., Fenton, W. H., Nelan, J. E., Patnaude, D. J., & Simon, F E. 

2001, ApJ, 562, 521 
Chaboyer, B., & Krauss, L. M. 2002, ApJ, 567, L45 
Clem, J. L. 2006, PhD thesis, Univ. of Victoria 
Clem, J. L., VandenBerg, D. A., & Stetson, P B. 2008, AJ, 135, 682 
Davis, D. S., Richer, H. B., Anderson, J., Brewer, J., Hurley, J., Kalirai, J. S., 

Rich, R. M., & Stetson, R B. 2008, AJ, 135, 2155 
De Angeli, F., Piotto, G., Cassisi, S., Busso, G., Recio-Blanco, A., Salaris, 

M., Aparicio, A., & Rosenberg, A. 2005, AJ, 130, 116 
Delahaye, F, & Pinsonneault, M. H. 2006, ApJ, 649, 529 
Demarque, P., Guenther, D. B., Li, L. H., Mazumdar, A., & Straka, C. W. 

2007, Ap&SS, 447 

Dotter, A., Chaboyer, B., Jevremovic, D., Kostov, V., Baron, E., & Ferguson, 

J. W. 2008, ApJS, 178, 89 

Dunkley, J., et al. 2008. larXiv:0803.0586l 



rial are those of the author and do not necessarily reflect the 
views of the National Science Foundation. 

Funding for the SDSS and SDSS-II has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and the 
Higher Education Funding Council for England. The SDSS 
Web Site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Con- 
sortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, As- 
trophysical Institute Potsdam, University of Basel, Univer- 
sity of Cambridge, Case Western Reserve University, Uni- 
versity of Chicago, Drexel University, Fermilab, the Institute 
for Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cos- 
mology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, 
the Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico State 
University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. 



Eggen, O. J., & Sandage, A. R. 1959, MNRAS, 119, 255 
Eisenstein, D. J., et al. 2006, ApJS, 167, 40 

ESA. 1997, The Hipparcos and Tycho Catalogues (ESA SP-1200) 
(Noordwijk: ESA) 

Fornal, B., Tucker, D. L., Smith, J. A., Allam, S. S., Rider, C. J., & Sung, H. 

2007, AJ, 133, 1409 
Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., & Schneider, 

D. P 1996, AJ, 111, 1748 
Geffert, M., & Maintz, G. 2000, A&AS, 144, 227 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 
Girardi, L., BertelH, G., Bressan, A., Chiosi, C, Groenewegen, M. A. T., 

Marigo, P, Salasnich, B., & Weiss, A. 2002, A&A, 391, 195 
Girardi, L., Grebel, E. K., Odenkirchen, M., & Chiosi, C. 2004, A&A, 422, 

205 

Gratton, R., Bragagha, A., Carretta, E., & Tosi, M. 2006, ApJ, 642, 462 
Gratton, R. G., Fusi Pecci, F., Carretta, E., Clementini, G., Corsi, C. E., & 

Lattanzi, M. 1997, ApJ, 491, 749 
Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 
Grundahl, F, Stetson, P B., & Andersen, M. I. 2002, A&A, 395, 481 
Gunn, J. E., et al. 1998, AJ, 1 16, 3040 
Gunn, J. E., et al. 2006, AJ, 131, 2332 

Gustafsson, B., Edvardsson, B., Eriksson, K., Graae Jorgensen, U., Nordlund, 

A., & Plez, B. 2008, ArXiv e-prints, 805, arXiv:0805.0554 
Hansen, B. M. S., et al. 2002, ApJ, 574, L155 
Harris, W. E. 1996, AJ, 112, 1487 

Hauschildt, P H., Allard, F, & Baron, E. 1999, ApJ, 512, 377 
Heasley, J. N., & Christian, C. A. 1986, ApJ, 307, 738 
Hogg, D. W., Finkbeiner, D. R, Schlegel, D. J., & Gunn, J. E. 2001, AJ, 122, 
2129 

Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 

Ivezic. Z., et al. 2004, Astronomische Nachrichten, 325, 583 

Ivezic, Z., et al. 2008a, ApJ, 684, 287 

Ivezic, Z., Tyson, J. A., Allsman, R., Andrew, J., Angel, R., & for the LSST 

Collaboration 2008b, arXiv:0805.2366 
Jensen, E., Boesgaard, A. M., & Deliyannis, C. R 2006, BAAS, 38, 1134 
Jimenez, R., & Padoan, P 1998, ApJ, 498, 704 
Johnson, H. L. 1957, ApJ, 126, 121 
Juric, M., et al. 2008, ApJ, 673, 864 
Kaiser, N., et al. 2002, Proc. SPIE, 4836, 154 

Kim, Y.-C, Demarque, R, Yi, S. K., & Alexander, D. R. 2002, ApJS, 143, 
499 

Kirby, E. N., Guhathakurta, R, & Sneden, C. 2008, ApJ, 682, 1217 

Korn, A. J., Grundahl, F, Richard, O., Barklem, P. S., Mashonkina, L., Collet, 

R., Piskunov, N., & Gustafsson, B. 2006, Nature, 442, 657 
Korn, A. J., Grundahl. R, Richard, O., Mashonkina, L., Barklem, P. S., Collet, 

R.. Gustafsson, B., & Piskunov, N. 2007, ApJ, 671, 402 
Kraft, R. P, & Ivans, 1. 1. 2003, PASP, 115, 143 (KI03) 



26 



An et al. 



Ki'aft, R. P., & Ivans, I. I. 2004, in Carnegie Observatories Astropliysics 
Ser 4: Origin and Evolution of tlie Elements, 2004, ed. A. 
McWilliam & M. Rauch (Pasadena: Carnegie Observatories), 
http://www.ociw.edu/ociw/symposia/series/symposium4/pro ceedings.htmli 
Ki-auss, L. M., & Chaboyer, b. 2U0j, Science, 2yy, bi 
Lejeune, Th., Cuisinier, E & Buser, R. 1997, A&AS, 125, 229 
Lejeune, Th., Cuisinier, E & Buser, R. 1998, A&AS, 130, 65 
Lupton, R. H., Gunn, J. E., & Szalay, A. S. 1999, AJ, 118, 1406 
Lupton, R. H., Ivezic, Z., Gunn, J. E., Knapp. G., Strauss, M. A., & Yasuda, 

N. 2002, Proc. SPIE, 4836, 350 
Michaud, G., Richard, O., Richer. J., & VandenBerg, D. A. 2004, ApJ, 606, 
452 

Mihalas, D., Hummer, D. G., Mihalas, B. W., & Daeppen, W. 1990, ApJ, 350, 
300 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Oiiglia, L., Valenti, E., Rich, R. M., & Eerraro, E. R. 2006, ApJ, 646, 499 
Ferryman, M. A. C, et al. 2001, A&A, 369, 339 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, S. M., Lupton, 

R. H., & Ivezic, Z. 2003, AJ, 125, 1559 
Pinsonneault, M. H., Terndrup, D. M., Hanson, R. B., & Stauffer, J. R. 2003, 

ApJ, 598, 588 

Pinsonneault, M. H., Terndrup, D. M., Hanson, R. B., & Stauffer, J. R. 2004, 

ApJ, 600, 946 
Reid, I. N. 1997, AJ, 114, 161 
Reid, N. 1998, AJ, 115, 204 

Renzini, A. 1991, NATO ASIC Proc. 348: Observational Tests of 

Cosmological Inflation, 131 
Renzini, A., et al. 1996, ApJ, 465, L23 

liider, C. J., Tucker, D. L., Smith, J. A., Stoughton, C, Allam, S. S., & 

Neilsen, E. H., Jr 2004, AJ, 127, 2210 
Rogers, E J., & Nayfonov, A. 2002, ApJ, 576, 1064 
Rosenberg, A., Saviane, I., Piotto, G., & Aparicio, A. 1999, AJ, 118, 2306 



Salai-is, M., Chieffi, A., & Straniero, O. 1993, ApJ, 414, 580 
Salai-is, M., & Weiss, A. 2002, A&A, 388, 492 
Sandage, A. 1970, ApJ, 162, 841 
Sandquist, E. L. 2004, MNRAS, 347, 101 

Saumon, D., Chabrier, G., & van Horn, H. M. 1995, ApJS, 99, 713 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 

Schlegel, D. J., Einkbeiner. D. P, & Davis, M. 1998, ApJ, 500, 525 
Sills, A., Pinsonneault, M. H., & Terndrup, D. M. 2000, ApJ, 534, 335 
Smith, J. A., et al. 2002, AJ, 123, 2121 

Sollima, A., Beccari, G., Eerraro, E. R., Eusi Pecci, E, & Sarajedini, A. 2007, 

MNRAS, 380, 781 
Spergel, D. N., et al. 2003, ApJS, 148, 175 
Stetson, R B. 1987, PASP, 99, 191 
Stetson, R-B. 1994, PASP, 106, 250 
Stetson, R B., & Hanis, W. E. 1988, AJ, 96, 909 
Stetson, R B., Vandenberg, D. A., & Bolte, M. 1996, PASP, 108, 560 
Straniero, O., Chieffi, A., & Limongi, M. 1997, ApJ, 490, 425 
Stoughton, C, et al. 2002, AJ, 123, 485 

Tucker, D. L., et al. 2006, Astronomische Nachrichten, 327, 821 
VandenBerg, D. A. 2000, ApJS, 129, 315 
VandenBerg, D. A., & Clem, J. L. 2003, AJ, 126, 778 

VandenBerg, D. A., Richard, O., Michaud, G., & Richer, J. 2002, ApJ, 571, 
487 

Vandenberg, D. A., Stetson, P B., & Bolte, M. 1996, ARA&A, 34, 461 
Venn, K. A., Irwin, M., Shetrone, M. D., Tout, C. A., Hill. V, & Tolstoy, E. 

2004, AJ, 128, 1177 
York, D. G., et al. 2000, AJ, 120, 1579 
Zinn, R., & West, M. J. 1984, ApJS, 55, 45 
ZoccaU, M., et al. 2001, ApJ, 553, 733 



